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Itv was during World War ‘that r radar is is an excellent nt instru- 
ment for obtaining dé data on the ‘Spatial distribution of precipitation. During 
the past ten’ years considerable research has been supported by the Depart- 
ment of Defense at various institutions to determine the application of radar 
ine Meteorology. _ Efforts have been c concentrated dont the quantitative measure-— 
"ment of precipitation and on the use of radar for the detection and tracking of 
“a tornadoes, severe winds, and hailstorms. As a result of studies atthe _ 
_ Massachusetts Institute of Technology, th the Illinois State Water | Survey, the 
U. S. Weather Bureau, and the Department of Defense, networks of radar are 
now being installed by various agencies interested in weather. In some cases 
‘radar | sets are being procured by industry to improve short- -range ‘een 
of precipitation and weather which are pertinent to their operations. am 
; a Illustrations show that radar occasionally detects the initiation of a new i 
q precipitation area one to six hours before it reaches a regular hourly r re- hor 
porting weather station. areal extent of precipitation is fairly accurately 
_ displayed and the relative intensity of precipitation can be determined. 
_ Changes in the intensity and nature of the precipitation as well as in exact a 
timing of the event are readily determined. Precise quantitative measure- 
poner for all types of rains have not been obtained to date. Despite this 


4 


radar to improve their skill in 1 short - ~range forecasting of rainfall 
river-stages especially during flash floods. 


By observing severe thunderstorms on radar scopes that are associated 
_ with tornadoes, the meteorologist is learning to track existing storms and _ 
perfect his short-range forecast of areas for future warning. _ The loss of life 


‘during the 1957 tornado season was oer reduced through use of radar to — a 


Note: Discussion open until June 1, 1959. extend on one month, a 
written request must be filed with the Executive Secretary, ASCE. 4 Paper 1901 is 
_ part of the copyrighted Journal of the Hydraulics Division, Proceedings of the Mn 
American of Civil Engineers, Vol. 85, Ho. HY 1, January, 1959. 


. Presented at the ‘meeting g of the Hydraulics Division o: of the ASCE in 1 Chicago f- 


"on February 25,1958. 
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INTRODUCTION 


all of areal rainfall and often hourly or storm data are not available con-— 
current with the event to river forecasters, meteorologists, and other ~ ; 
a rainfall data. The raingage merely measures the rainfall in the immediate 7 
_ Vicinity of the collector. Numerous investigators(1,2) have shown the inade-_ ane 
‘quacy of point measurements. Figure 1 is an example from one of the Illinois 
State Water Survey’ s raingage networks which illustrates the variability « of 
storm, monthly, and growing season rainfall for a 400-square mile area in 
> al ~ central Illinois. Records from fifty recording raingages were used in — "a 
ing the detailed isohyetal patterns. ¥ 
_ During World War II it was observed by 1 military personnel that certain 
of radar are detectors of precipitation. Preliminary studies 
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‘Figure Daily, monthly and seasonal | rainfall for 
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"RAINFALL MEASUREMENTS 


‘In 1948 the Illinois’ State Water Survey became in the utility 
radar for quantitative measurement of summertime precipitation in Illinois. | 

; Networks of closely spaced raingages were established to define the repre- 4 

sentativeness of rainfall observations and compare | radar - 
with gage-measured point and areal rainfall.(3) 

Numerous studies have been n completed in the eight years to Ping 

4 mine the utility of radar in measuring rainfall. Preliminary findings of the 
various research organizations are presented in this paper in order | to (o sum-_ 

marize progress to date. _ Difficulties have been experienced because re- » - 
search on quantitative precipitation measurement by radar have had to >a 


q satisfactory test equipment. Despite several limitations, “many ‘applications 


of f radar are 8 


: energy which. is focused into a narrow beam of invisible energy by a rotating 
antenna, much like a search light. _ This pulse of energy travels at the speed 
of light. If the beam strikes an object, such as an airplane, a building, ora 
group of raindrops in a rain cloud, a small portion of the energy is reflected _ 

: back as an “echo” to the point of transmission. The return signal is ampli a 
fied and presented on a cathode ray tube, which is referred toas a ‘Plan- or 
Position-Indicator. The range and bearing of the object are automatically > 

~~ ‘Research it in 1 this field by previous investigators(4, 5) indicated that radar 

echoes from rain clouds are the result of back scattering - of radio ) energy by — 

4 raindrops falling through the atmosphere. The theoretical expression that _ 

a has S been developed indicates ‘om the energy (power) received from raindrops 


: where Py i is received power, Py is power transmitted, r is the di stance > 
trom the radar set to the reflecting raindrops, N is | the number of — 
per unit volume, _d is the diameter of the raindrops, and K absorbs a con- 
oe for the refractive red of w water and ; and a number of constants w which a are 


parameters of the radar set. yy a Bite the 

does not measure rainfall but measures the re- 
flectivity, Z, which is pte to 2 Nd®°. Several expressions for rain- 
fall rate in terms of 2} Nd6 have be raindrop size data. data. 


“When one of these expressions, 

logarithmic form is, for ‘the APS-15 radar eystem, 
3 


is 
7 
a: 
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_ 4 Before considering any of the 1e research t to date, it may be instructive to 
consider - some of the assumptions that are ‘made. _ The theory | is based on a 
scattering by spherical raindrops whose diameters are much smaller than an 
the wavelength of the radar. The absorption of the signal, also known as “- 
~ tenuation, and the back-scattered energy a are both increased by deviations ‘bd 
- from non-spherical drops. Most drops larger than 1.8 mm are non- spherical 3 
_ Another basic © assumption | is that the radar beam is 3 filled with ne 


—, is not always correct for storms at long ranges, since the rain- 
storm may be ‘smaller than the cross section of the beam. | Also, since the 
height of the beam under normal refraction conditions is proportional to the 
_ square of the range, it is evident that at longer ranges the beam is higher in : 
_ elevation and i therefore less likely | to be filled. At a range of 100 miles a 
height of the center. of the beam is approximately 5000 feet. Owing tothis 
limitation, present- -day radar used for quantitative rainfall research should 
usually be restricted to 100 miles. 

4 Another consideration of importance in quantitative precipitation measure-_ 
‘ments is that of the relationship of Z, radar reflectivity, to R, rainfall rate. 
Based. 1211 minutes of data, collected by the Illinois State Water Survey, 

using a a raindrop camera,(6) Table 1 shows how R varies for giv 

Z under different conditions of of drop 


Ze for Various 
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drops in camera, Jones(6) distance of 30 miles 


‘The fact that the is a some aloft while the 
precipitation and drop size measurements are determined at the surface is . 
problem of indeterminate | seriousness. . Due to their smaller terminal veloci- 


ty, ‘small drops are carried farther horizontally y than larger drops. This — 
ari may become serious ina wind field where the vertical wind shear is 


may occur under vertical wind shear conditions. Evaporation. and 
a of the drops may also change the reflectivity and rainfall rate 


between the radar volume and the ground. Thus, an error in surface rain- 
fall estimation can be induced by those physical processes. __ Nera ES 7 
_ A further consideration is stability of the radar. - Normally the receiver . 
Sensitivity i is S subject { to fluctuation and should be monitored quite frequently. 
ff Another rather serious limitation at 3 cm and shorter wavelengths is = 
: _tenuation; that is, loss of signal a due to intervening rain which ab- 
~ King(7) determined experimentally that ‘the . attenuation is proportional to the 
_ product of rainfall intensity and distance into the storm. McGill University(8) ui 
and the Illinois State Water Survey using observed drop size data and Mie 
= scattering theory found that attenuation per unit distance was related expon- 
4 entially to the rainfall rate. Robertson and King’s relation gives higher values" 
attenuation for rainfall rates below 90 Above 90 ex- 


_ As ‘Part | of a 2 project at the Illinois State Water Survey, ‘sponsored by can 


“fulness of radar in estimating quantitative precipitation, areal rainfall _ 
for one ‘minute were e compared with areal mean radar 0 one- minute 


“gages in a 100- -square mile area located 20 miles northwest of the radar sta- 
tion. Radar- indicated amounts were obtained from a calibrated APS-15A, a 


radar Areal means were determined by planimetering the radar 
‘Figure 2 shows the comparisons between radar power ratio and ee 
amounts. The data(9) were derived from 95 minutes of summertime rains — 

where no echo intervened between the radar site and the network. Cases of — 
marked dissimilarity between 1 radar and raingage patterns were also elimi- 
nated from the analysis. A lag up to 2 minutes between radar and raingage 
amounts was allowed to permit the precipitation to fall to the g: ground from a 
beam height. Considerable scatter inthe datais evident. 
lan From recent analysis of an unusually heavy rainstorm, (10) the sia 

conclusions regarding quantitative re radar ¢ observations at 3 cm were re pre- 


“The CP5- 9 radar est: is satisfactory for portraying the areal extent of 


i _ edge and the horizontal extent of an approaching storm in the absence of in- “4g 
_ tervening echoes, as shown in Figure 3 for t the October 9, 1954 | storm ov over 
northern Ilinois. At 1730 CST numerous rain cells are evidenced in areas 


between recording raingages. _ However, at 20 2030 30 CST neither the depth 
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Figure 2. Relationship between rainfall int intensity 
and radar- -received power factor. 
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nor the intensity of a storm is accurately portrayed. Consequently, objective 
methods of obtaining accurate quantitative rainfall estimates with 3- cm radar 
do not appear feasible. However, an experienced radar meteorologist, using 


synoptic weather information and climatological relationships, along with 


_ knowledge of the limitation of 3-cm wavelength radar, is able to make better : 


estimates of the rainfall in an area than can be obtained from analysis of nod 


Synoptic weather maps alone. he’ 


an In another study the r. radar change in n intensity | was s compared with the rain- a 
«gage intensity change. Using 5-minute amounts for cases of insignificant at- _ 


Robertson- -King’s equation, the Illinois State Water has 
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‘ama 3. Radar echoes and hourly rainfall amounts for October 9 ae 
1954 over northern Illinois and Indiana. (solid 
full gain; solid areas, 26-DB reduction. All station 
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"January, 1959 
strikingly dissimilar patterns were eliminated. Results 
based ona comparison of point rainfall indicate no discernible improvement mm 
in the 3-cm radar’s ability to measure rainfall. Areal mean rainfall compar- 
isons were | made | in an effort t to reduce | observational error, but r no Sig - 


efforts with attenuation correction indicate that other factors are important 


attenuation. Some of these discussed earlier 
y; 


_ The Weather Radar group at the Institute of 
_ under sponsorship of the U. Ss. Army Signal Research and Development v 
Laboratory, has been testing its pulse integrator with a 10-cm radar as >» 
method of determining intensity of rainfall. The pulse integrator is an elec- 
_ tronic signal intensity averaging device. Utilizing the ranging unit of t the ‘ 


radar, the peak voltage of the return signal at a point or along an azimuth is 


electronically averaged. Gating allows the scope to be divided into numerous _ ¢ 


from a MIT report(11) is an example of the type of information readily availa- 
ble. Note the similarity in the data when a radar a areal Te is anne 


areas so that rainfall values can be determined. Figure 4 which was prepared . 


with a point raingage rate. 
_ Correspondence with MIT personnel indicates that satisfactory aati | 
(approximately 50 percent) is obtainable, provided the relation between signal 
intensity on their 10-cm radar and rainfall rate is obtained from a series of 
7 measurements rather than from purely theoretical considerations. Also, one 
or two sensitive raingages in the region being covered should be used as -* 
check in any storm. Figure 5 shows the results of some recent work at ‘- 
MIT.(12) For light rainfall rates of less than 1/2 inch per hour, there isa | 
similarity in traces between the radar and a raingage » Ho However, an error of 


_ Texas A & M in cooperation with the U. S. Weather Bureau has 0% inter- 


point out the well-known fact that a storm moving p eet sien a basin ; 
produces a much sharper hydrograph than one moving up the same basin, al- 

_ though the actual volume of rain may be the same. They recommend deter-_ 
mining the “velocity of maximum concentration v vector” for the individual 
i - drainage basins. + This is defined as the direction and speed of movement of i 7 

_ given storm that will produce the maximum crest, other things being equal. " 

& ‘Radar scope composite photos a are made by photographing the | scope at 5- — 
= intervals for an hour and ‘making one print from these 12 + Rha 


% — isohyets for the 3-hour period were drawn to fit these data. ‘Plone 6 
shows the results of this analysis. The assignment of values to these radar — 
ie isohyets, based on the 3-hour recorded precipitation amounts, was attempted 
but proved difficult and in many cases grains radar isohyet would ee 
too wide 
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q hourly raingage data. Radar-indicated rainfall amount 


PULSE INTEGRATOR (mm/h) = 6 


Comparative pt pulse integrator or and rainfall values for 


July 10, 1952 (courtesy of Massachusetts 


» radar isohyets are labeled as A, B, and C. At tir At core 


x The University of Miami in cooperation with the U. S. Weather Bureau is | 

comparing 2-hour composite photos with observed rainfall values, using an 

SP-IM, a 10-cm radar set. At 10-cm there is the advantage of less attenua- 
tion than at 3-cm. Since the radar scope, also known as the Plan- cn 
Indicator (PPI), is an intensity modulated scope, the more intense areas over a 


the period of integration will produce greater darkening of the film oe anes 
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nse ar . Also, a storm of longer duration will expose ine: film more 
than one of short duration. _ After many types had been tested, Kodak- © 
Ektalure film was selected owing to its slow speed and wide latitude of - 
1 Many photographic and radar problems are involved in measuring g rainfall. 
_ quantitatively by radar. Critical control of the photographic developing __ 
process, film shrinkage, and halation are photographic problems. At close | 
_ range, the scope becomes saturated. A logarithmic receiver instead of the 7 
conventional linear receiver is one solution to this problem. 7 The fact that 
Z the velocity of the sweep in azimuth is proportional to range causes echoes at 
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close range to appear brighter than those farther away. This effect is being — . 
_ 3 compensated for by a filter, the density of which is inversely proportional to = : 
a range. Furthermore, critical control of radar performance and setting must — 


re 7 shows an opacity-range- -rainfall nomogram prepared from these 
a 15) Although the scatter on these diagrams is large, a potentially use- _ ¥ 
ful relation is indicated. This is a rapid and inexpensive a of ee 
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a Future developments will probably inc include cos photographic i inte- 
: gration. . This is a relatively inexpensive and quick way to obtain areal mean 7 
Fainfall over several hours. More use is likely to be made of color presenta- 
tion. Different colors could be used to represent a variety of intensities and, 
thus, enable an operator to judge the intensity distribution at a glance. _ Stor- — 
age and memory tubes could enable the operator to survey the storm in time 4 
n and space. Electronic computers may also be used but they are expensive. a 

_ Transmission of an actual scope picture to distant places may ultimately i 
replace the coded message of today. Work is already in progress on trans- - 
mitting the scope picture by facsimile, employing a low frequency to permit _ 
use of land ‘Hines. Th Thus, large scale composites can be made at an analysis 
a the radar sets themselves we can cabbih the development ofa greater 
dynamic range, , which will eliminate scope saturation from the ———— : 
=a integration process. q Circuitry that will compensate for range attenuation — 

may be developed. Constant height PPI presentation may add to a better < 

understanding of the precipitation mechanism. = = 


4g Development of equipment that would make radar information more ac- 
-_cessible, useful, and reliable is in prospect. Analysis will most likely be 
aimed at determining areal precipitation over a several hour period and 
developing models of patterns. 


Severe st ‘storm Tracking by Rada Radar » - 


- Since 1953 radar has proven itself to be a valuable instrument in detecting 
‘the a areas of severe hailstorms and tornadoes. Incidental to Illinois 


plete radar recording of the development of a tornado causing "damage esti- _ 
mated at $3,000,000 was m made. (16) Figure 8 illustrates the | PPI for April . 
1953, “showing a hooked appendage ‘developed on the southwest corner of “a 
_ large thunderstorm. A tornado was centered under the vortex part of the 
‘ appendage. Since then, numerous observations have been recorded by radars ° 
te in Massachusetts, Oklahoma, ‘Texas and Kansas. It appears that the maximum — 
detectable range of hooked-shaped echoes is about 50 miles. Actually, the a 
parent thunderstorm cloud may be detected | at ranges up t to 250 miles. RAS = . 
@ Exact data on the value of radar through better warnings during local A 
severe storms is difficult to ascertain, since data on past storms are very” 
"incomplete. _ Experience has indicated that radar is very effective in in the de- a 
tection and tracking of severe storms. ‘The network of radar stations needs — 
_ to be expanded and improved upon. Better communications between the Oob- 
Server of the radar scope, the forecaster, , and the public should be estab-— e 


lished. Eventually, the loss of life from tornadoes will be negligible. 


that point or areal measurement of precipitation by radar are feasible, | 
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_ Figure 8.. Photos of radar PPI with 10-mile range markers = 
showing development of cloud pattern associated _ 


swith the tornado of April 9, 1953, 


RAINFALL MEASUREMENTS 
to discrepancies of considerable magnitude at accuracy 
_ seems to be limited by the uncertainty in the empirical relation between sig- 
nal intensity and rainfall rate which are based upon several assumptions. — 
os , one must keep in mind that rainfall amounts from raingages are 
often not representative of areal mean rainfall. 
The applications of radar in hydrometeorology are numerous. Concurrent 
with the event, the detection, areal extent, relative intensity and changes in 7 
; intensity, duration, ‘storm velocity, and orientation of precipitation with /~ 


respect to the basin are readily recorded for the immediate use of the river : 


or the synoptic forecaster. 


a Considerably more research and development of equipment should be pur- 
sued to solve some of the problems currently confronting the radar meteorolo- 


4 Radar data will become more valuable as a network of stations is b-pwering 
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To meet ‘the increased needs of ocean navigation, ‘the Federal project for 7 
improvement at the mouth of the Columbia River was modified by the River — 
12 Harbor Act of September 3, 1954, to provide for an entrance channel a 
‘suitable alignment having a depth of 48 feet at mean lower low \ water and a 
- width of one-half mile. Under the modified project, this channel would be ob- 
|Z by dredging and the construction of a spur jetty on the north shore, in 
_combination with existing jetties and contraction works, at an additional cost 
_ currently estimated at $10, 340,000 for new work. The difficulties associated — 
— with prosecution of this modified project call for a re-evaluation of the TT 


al forces controlling the regimen of the entrance area. tt 


4 _ The Corps of Engineers is attempting to secure by vigorous dredging a 
‘channel at the entrance to the Columbia River with least depth of 48 feet and 
of such width as may be permissible within the limits of available funds and 
- dredging plant and the minimum requirements of present navigation. Although: 
ee depths of only 35 feet are authorized and maintained along the main _ 
_ Columbia River channel from the entrance to the mouth of Willamette River — 
and thence to Portland, a distance of 117.5 miles, an additional usable depth _ 
- of 13 feet through the entrance is necessary to assure safety to ocean vessels 
_ particularly during the heavy sea swells and storms which occur at that lo- a 14 
a cality. With this full depth assured, the physical capacity of the entrance to — 


carry ocean commerce will be practically unlimited. However, maintenance ~ ; 


. Note: Discussion open until June 1, 1959 . To extend the closing date one month,a 
a written request must be filed with the Executive Secretary, ASCE. Paper 1902 m 
- part of the copyrighted Journal of the Hydraulics Division, Proceedings of ood Me 

American of 85, No. HY 1, 1959. 
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project and the work of securing the new project di- 
- mensions have been made difficult, not only bya adverse weather and sea “i 
ditions, but also over the years by the growth of the hazardous Clatsop _— 
‘shoals. These shoals, lying along the channel side of the South Je Jetty, have 
- gradually forced the navigation channel toward the North Jetty, as ‘shown 4 

‘F igs. 3 and 4, inc reasing annual dredging requirements and creating a1 “sg te 
these difficulties, al usable. channel, although less than full di- 
———_ has been maintained. _ During the period 1946 to 1956, ocean traffic ; 


_ million tons. per year, as s shown | in 1 the following tabulati tabulation: 


MILLIONS OF TONS 


of a entrance with depth of 4-1/2 fathoms between 
Disappointment and Point Adams. In 1839, Sir Edward Belcher revealed aad 
7 existence of two channels at the entrance separated by a large shoal area — 
which forced the deeper channel northward against | Cape Disappointment to. 
turn 2 abruptly southward and join a generally shallower channel along Point 
7 Adams. _ This shoal area, called the Middle Sands, grew in size until its 18-_ 4 
foot depth « contour, in 1874, embraced a very large a area seaward and south of | 
Cape Disappointment, and covered almost all of the area now occupied by the 
present entrance. A portion of the Middle Sands formed an island inside the 
entrance which later separated f from the main shoal area and, ‘migrating to. 
the north to form Sand Island, restored a single entrance channel to the mouth 
of the Columbia River in 188 5. At that time, _the main portion of the Middle 
Sands had moved westward | to form a southward extension of Cape Disap- a 
eo and ma depths of 20 feet prevailed on the ocean bar as shown _ 
a es to recommendations of the Board of Engineers for the Permanent ij 
_ Improvement of the Mouth of the Columbia River, (1) Congress on August 2, > 1 
1882, adopted the original permanent navigation project for the Columbia __ - 
aD River entrance. This project provided for a bar channel 30 feet deep to be ob- _— 
tained by construction of a South Jetty, 4-1/2 miles in length, lying in a gener- “4 
ally northwesterly direction from its base at Point Adams. Construction of ‘ 
the jetty was started in April 1885 and was completed. with four groins along 
its north side, in October 1895. At that time the crest of the structure sloped 
% from 12 feet above mean lower low water at its base to plus 10 feet at a a 


a t -1/8 miles from the shore, tence to plus 4 feet : at its outer end. A total of 
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946, 000 tons of rock was s placed in the original construction of me South Jetty 

‘ Although shoal 4 areas I marking the birth of Clatsop Spit awe to form along © 

the north side of the structure during its construction, a general imgrovement ‘- 

of navigable depths over the bar was evident in this period. By 1895 mini- — i 
mum depths of 31 feet were available over the bar on an which, 
_ under the influence of the jetty, had migrated to the north about three miles. — 


migration of the channel to the north; by 1896 near closure of the 30- -foot 

depth contour on the bar occurred; by 1898 depths of only 29 feet were avail- 
able; by 1901 entrance depths had deteriorated further; and in 1902, instead of 
one, there were three separate entrance channels, each being about 22 feet _ 

by deep, as shown in Fig. 2. In the meantime, ‘Clatsop Spit ‘shoals had formed a 
hook-shaped island through the South Jetty at its mid- -length which later ex- | 
tended to form a solid land connection with Point Adams. 


March 3, 1905, adopted of the Board of on 
Project for Improv ement of Mouth of Columbia River. (2) These recommen-— 
dations provided for modification of the existing project to secure an entrance 
channel one-half mile wide and 40 feet deep by extension of the South Jetty, _ 
construction of a North Jetty and groins, and dredging. - Extension of the South 
Jetty, completed to a total length of 6.6 miles in August 1913, restored .. ; 


single entrance along its former alignment to the mouth of the river. By 1911 


a depth of 24 feet was again available across the bar. As completed, the South — 


Jetty had an average top width of 25 feet and an elevation above mean sea 
level of 10 feet from the shore to the “knuckle” end of the original project, 
from whence the height was increased by stages to an elevation of plus 24 ~~ 
feet at its outer end. Construction of the North Jetty began in 1913 and by the : 
next year, when only about 1 mile in length, entrance depths had inc reased to 
30 feet. As constructed, the North Jetty had a top width of 25 feet, side slopes ~ 
of 1 on 1-1/2, and crest elevation of 28 to 32 feet above mean lower low a 
Completion of the jetty to its full length of 2.4 miles in August 1917 was ac- 
companied | by further improvement of entrance depths to 40 feet, which depths 
were available through a narrow bar channel. By 1920, the 40-foot depth 
contours had retreated considerably providing an entrance channel about one 
mile in width; and , during | the next five years, these contours retreated farther 
providing a minimum entrance depth of 45 feet in 1925. The next two years 
have been recognized as the period when the most satisfactory channel con- 

4 ditions prevailed at the mouth of the Columbia River ‘By 1927, minimum = 


: - depths of 47 feet were available over the bar, as shown in Fig. 3. Itis a 
_ esting to note that for the first time, depths in excess 0 of 40 feet prevalies 
along the channel side of the North Jetty, 
_ However, concurrently with improvement of navigation conditions, ode 
_ developed behind the North Jetty and formed a solid land mass extending al- _ 
most to its full length. _ Clatsop Spit continued to g grow y until by 1931 its il ; 


westward and northward extension had reduced entrance depths to about . oo 

feet. Rehabilitation of the South Jetty was started in 1931 and completed in 

_ February 1936; then, later in 1936, the outer 500 feet of the structure was _ 
"impregnated with a hot asphaltic mix in an attempt to prevent raveling at this 

= by heavy seas. This expedient was not fully successful and a solid 

Concrete was constructed at the outer end in 1941, Further growth 
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be a These favorable conditions, however, were short-lived because of continued © 
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- deeper water along the channel side of the outer end of that structure. In an 
“effort to prevent continued movement of the channel to the north, four perme- 
able dikes were built along the south shore of Sand Island, and Jetty “A”, ex- 
tending southward from Cape Disappointment, was completed in 1939. The _ 

continued westward movement of Clatsop Spit resulted in the removal o of more re 
of the shoal area protecting the North Jetty thereby causing some damage . 
‘from undercutting along the channel face of the structure. During 1938 and — 
1939 the jetty was rehabilitated and a concrete terminal block was placed at 
its outer end. By 1944, although the 40- foot | depth contour - opposite Jetty <A” 
had moved southward, this local improvement in conditions was accompanied _ 
by marked growth of Clatsop Spit to the west and north opposite the end of the 
_ North Jetty. Continued advance of the spit through 1951 coincided with definite 
7 movement of the 40-foot depth contour along the north side of the entrance 
_ channel to the east and south and in 1952, for the first time, an inner bar in © 


_ extension of Clatsop Spit had formed between the outer end of the jetties. | an 
_ In addition to 410,000 tons of stone used for rebuilding the outer end ofthe _ 

original South Jetty, a total of 4,427,000 tons of stone was placed in the ex- 
tension of the South Jetty, completed in 1913, at a total cost of $5,657,000. 
Nearly 3 million tons of stone were used in the construction of the North Jetty 
a cost of $4,319 ,000. Reconstruction of the South Jetty involved the 
- placement of 2,288,000 tons of stone in the superstructure at a cost of ‘>: 
$4,930,000. | Rehabilitation of the North Jetty required the placement of 
244,000 tons of stone and 234,000 tons were used ene construction of ‘Jetty = 
“ar at a total cost of $1,271, 000 
_ As ample depths prevailed over the ocean bar, po under the 40- -foot 
entrance project was confined to the face of Clatsop eee dredged — 


since 1939 project were as follows: | 


FISCAL YEAR CUBIC YARDS YEAR CUBIC YARDS 


1,319,800 


The full project width was obtained d only during the last three fiscal years 
on April 15, 1956 initiated the work of dredging the 48- “foot entrance channel; 
and, by the ‘end of the « dredging ‘season in October, minimum depths of 48 feet. 
Ss the full channel width of one-half mile through the inner and outer bars 
were obtained. These operations were accompanied by varied amounts of 
= apparent scour and shoaling. Heavy shoaling of the channel during the a 
dredging period between October and April, together with the lack of sufficient 
dredging plant, resulted in the decision to concentrate the 1957 dredging © —t 


_ program to the securing of project depth throughout a channel of 1500-foot _ 
= along the adopted alignment. These later operations were also ac- 
‘companied by natural | scour and pts and i heavy shoaling again occurred = 


of Spit ‘forced the navigation yn channel to the North Jetty q 


1939 

= 

— 


wn in Fig. 4, ‘The following table shows the results of the 1956 and S, 


1957 dredging programs and changes occurring within the adopted channel 
since initiation the as revealed by periodic condition — - 

4 DREDGING PROGRESS, MOUTH OF COLUMBIA RIVER, "FOOT FROJECT 


QUANTITY OF QUANTITY OF APPARENT SCOUR OR 
IN MATERIAL —DURING PERIOD, 1000 q 
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‘on During initial construction of the South Jetty, shoal areas, which a 
_ formed a land connection with Point Adams, gradually developed along the = 

outside of the original structure. Almost concurrently, shoaling occurred _ 

_ren the inside of the jetty forcing the main channel to the north. This a 

shoaling intensified with the extension of the South Jetty and with construction — 

of the North Jetty so that a large land mass, Clatsop Spit, now lies inside and» 

_ along the South Jetty. This land mass and its surrounding shallow depths — 

present a definite hazard to navigation, During the period January | 1947 sail 

- January 1952, the U. S. Coast Guard has reported a total of five vessels ol 
“* grounding in the entrance. In addition, owners of tankships carrying petrole- - 

_ um products reported that two of their vessels had cracked in crossing the + 
entrance and presumably had touched bottom. Although depths greater than 
48 feet are available ina natural channel inside the outer bar, just northwest : 

of the Clatsop Spit shoals, as shown in Fig. 4, this channel is particularly _ 
hazardous during severe storms and has a greater deflection angle than the | at 
- adopted channel that encroaches on the shallower depths of the spit. It has 
been possible to maintain project depths throughout usable widths : during a 
: portion of the time along the adopted entrance alignment by vigilant dredging 
_ when weather and sea conditions permitted, but such work has been costly. ; 
The problem resolves itself to one of securing the —_ project di- iv te 


Relating t to othe Problem 
‘The problem o of securing g dependable depths along a suitable alignment at - 
“tne mouth of Columbia River has presented a challenge to the Corps of Engi-_ 
neers for many years. complexities associated with this problem were 
=. _Fecognized as early as 1882, when the minority report of the Board of Engi- - 
‘neers for the Permanent Improvement of the Mouth of the Columbia River(1) 
and alternating currents, and the prevalence of sand in 1 the vicinity, in ape 
_ positions, introduce obstacles of so serious a nature that they may be said to 7 
_ place limitations upon the applications of principles at the Columbia River | 
Later, in 1903, the Board of on Project for of 


time to have a qt on the achievement of the deathet proj ect depths. — 
Subsequently, in 1932, an extensive program of current measurements across 
one principal river range, about 5 miles above the outer end of the jetties, in 
combination with other measurements including bottom sampling was omg 


reportt® on the » mechanical and petrographic characteristics of material 
. forming the river bottom, estuary, and beach at the mouth of the Columbia 
River was ‘prepared by Dr. Edwin ‘Hodge, Consulting Geologist, in April 
1934, ‘These : reports; studies under taken by Dr. M. P. O’Brien; -(5) the Tidal 
Model Laboratory at California, (6) salinity observations by the 
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Corps of phenomena observed by other Federal and State agencies, the 
for present knowledge of forces « controlling the regimen of the entrance area. 
From these data there has been extracted the following information ion pertaining 7 


to factors which may have a a significant bearing on the problem: fae 
The pattern c of lower Columbia River is influenced by 
- spring snowmelt runoff from upland areas and, to a lesser degree, by rainfall 
_ during the winter months. The low-flow period generally occurs during the | 
_ months of September through November and during this period river stages or 
may infrequently be as low as zero, the adopted low water stage at Vancouver. 
A this location, winter rains from December through February usually result | 
in gage readings | of 4 4 to 5 fi feet. Beginning in March, melting | snows gradually — 
: increase river stages to greater than 20 feet at Vancouver during ‘the latter > 
is _ part of May and through the month of June. River flows then gradually de- 
cline | to the low- flow period of September through November. Near ‘the mouth, 
= extreme low river daily flow is estimated at 60,000 cfs, the average flow a 
at about 255,000 cfs, the average spring freshet discharge at about 625,000cfs, 
and the maximum discharge at about | 1,300,000 cfs. Occasionally, winter 
4 freshets along the lower reach of the river exceed the average spring freshet. 
7 The average low water slope in the estuary is about 0.05 5 foot per mile. 


Tides” 
[ Tides at the mouth of the Columbia River have the diurnal inequality typi- — 
— cal of the Pacific Coast with a long ‘runout from higher high to lower low 
water. _ The mean range of tide is 6.5 feet, the range from mean lower low 
water to mean higher high water is is 8.5 feet, and extreme tides range from > 
‘ minus 2.6 feet to plus 11.6 feet mean lower low water. . Tidal effect extends 
_ upriver to Bonneville Dam, 140 miles, during periods of low flow. The tidal Be 


_ prism has been estimated at 782 square-r -mile feet. _ The current survey of of 
the he following characteristics of tides at the entrance: — 


‘TIDAL CHARACTERISTICS, OF COLUMBIA RIVER 


AVERAGE RIVER AVERAGE DURATION OF 
“MONTH DISCHARGE __TIDES IN HOURS _ 
1932000 CRS FLOOD EBB‘ FLOOD 


6.90 1, 
5.43 6.90 1, 
55 


“usually extends from October to and storm winds, generally 
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in the south quadrant, gradually veer to the ‘southwest and finally blow them a 
elves out as they move around to the west or northwest. Storms ar are often | the 
 # continuous for several days and are accompanied by extremely heavy seas : ‘sti 
mostly from the southwest. ‘During the summer season, the conditions se 
and sea are much less severe than those of the winter season. 
May through September, winds from the north or northeast begin in 
morning, increase in force during the day, and diminish at night. Occasional- — ort 


ly, during the winter and spring, there are several days when an ied -_ 

predominates and the weather is cold. 


Waves in excess of 40 feet in height have been reported by the Coast G G juard fe 

at the mouth of the Columbia River. Observations of wave conditions between en 

une 1934 and May 1935, as reported by Dr. M. P. O’Brien, indicate that the fo 

wave direction is the st. Weighted percentage of 


DIRECTION 


Southwest BA 


Northwest 


“Although h the C River is not ee a -sediment- -bearing 
stream and is relatively clear for the greater portion of the year, it is known _ 
that some material is carried in suspension during freshets. Measurements - 
below the mouth of Willamette ‘River made in sired and a at ‘Vancouver in 1 1950 


- the average freshet. This, however, does not include the bed load ofc coarser 
material which, available evidence seems to indicate, travels in waves that _ 
range from 8 to 14 feet high near the mouth of Willamette River but decrease © 
in height downstream and have crests from 300 to 800 feet apart. Tests on 


veal that these : waves ‘may move e from 30 to 50 fe feet ina ot in a single ales 


‘The significance « of the littoral to the problem may 
be evidenced by the tremendous build-up of material behind the North m may 
; which began with initiation of construction in 1913 and kept pace with the jetty 
as it was pushed seaward. . About 48 million cubic yards of material have 
been deposited in that area and, except for the part along the beach below high | 
tide, is now stabilized by a natural growth of grass and weeds. The fact that 
there has been no o similar build-up of materials behind the outer portions — 
the South Jetty appears to point to a predominate littoral drift to the south ~ 
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"undertaken in 1932 indicated that mean ebb tide velocities were iinian than 
- those of the flood tides for all stages of the upper river, , and | during freshet 7 
stages were much stronger. highest velocity ‘observed ¥ was 12.28 feet per’ 
second on the ebb tide at one- -tenth depth. These observations appeared to _ 
indicate that at high river stages, ebb bottom velocities were much higher — a 
than flood bottom velocities for all percentages of time. At intermediate _ 
river ‘Stage e, bottom veloc ities were lower than at high a with ebb — 
stage, ‘flood bottom velocities were generally: Considering the 
magnitude of ebb bottom velocities as compared to flood bottom velocities a ye 
and | their ‘greater transporting power, , it appeared evident at that time that the 


ebb flow was sufficient to maintain a channel between the 1e jetties of at least 40 


uard 4 feet deep at mean lower low water and that no further contraction of the _ : 
veen _ entrance appeared necessary for satisfactory channel maintenance. The 
following is a summary mean velocities measured in 1932 opposite Clatsop 
pit at strength oftides: 


Ebbs (Seven tides) 


Floods (Seven tides) 
seage 


Ebbs (Eight tides) 

Floods (Eight tides) 


‘At High Stage (May) 7 

nts Ebbs (Seventeen tides) 

— undertook a series of observations at the entrance use of the 
oi Radio Current waged to measure current velocities and direction at depths _ 
——e near the surface, at mid-depth, and near the bottom. _ Measurements were — 

on - | - taken at a total of seven stations, four located along the entrance channel at 

re- least one mile outside the jetties, one about two miles south of the outer =~ 
nt ‘, of the South Jetty, and two just inside ' the jetties. Although these measure-_ 
ments were made under difficulties and, accordingly, are incomplete, it = 
Me Re yeeren. to note that, while measurements at depths near the surface and 
may _at mid-depth indicated a general movement of water into the estuary on ann 
tty tides and out of the estuary on ebb tides, measurements at three and possibly \ 
jetty ‘tae stations appeared to indicate at the lower levels a general a of 7 7 
4 water into the estuary evenonebbtides. 
high In January 1958, further observations of current direction and velocity, 
that if salinity, and temperature were made by the U. S. Navy Hydrograph Office in 

s of 7 te -and adjacent to the Columbia River entrance _ Measurements made at a point 
about four miles upstream from the outer ends of the jetties during a river = 
“Eel discharge of 175,000 cfs indicated strong density currents. Upstream veloci- 
ee ties at the bottom ranging from 0.6 to 4.4 feet per s second were observed oeer- 
times of downstream and slack currents at the surface. 
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Measurements by the Corps of Engineers(7) in 1936 of the specific gravity | 
of surface, , mid-depth, and bottom samples at at 21 stations along the lower 22 7 
miles of the Columbia River e estuary indicated a definite difference in salinity 
between surface and bottom waters. _ During periods of flood tide on April 9 | 
and 10, when the river flow \ was about 150, 000 cis maximum observed 


1.024 along the bottom. “During periods of half to almost full tide on August 5 
and 6, when n the r river flow was about 165 000 cfs, the maximum 1 observed — 3 = 
densities were 1. 007, 1.020, and 1. .023, respectively. ‘These observations a 


= 


a: peared to indicate that on these dates the upstream limits of salinity pene- | 
trations were e about 15 15 and 11 miles above the outer ends of the jetties, re re- 
_ A number of salinity measurements were made by the ‘Oregon State College 
and the Oregon Fish Commission(8) on July 25, 1955 along the navigation La 


channel in lower estuary area 2 with findings as shown in 1 the following tabu- 


(MILES ABOVE MOUTH) 0/00 | PARTS. THOUSAND) AND DEPTH 


— 


1 (2.6 MI). 


‘ 

29.0 


00 


Lower high tide . occurred a at 0528, higher low. tide at 1138, and nidbee high 
tide at 1755 hours on 234 River: 000 cfs on that date 


1 
OBSERVATIONS, COLUMBIA RIVER ESTUARY, 1955 
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= the 1 following mies these agencies conducted salinity o observations at one 
station in mid- channel off Point Adams Coast Guard Station Station, ‘Mile’ 7. 5, 


"SALINITY OBSERVATIONS NEAR POINT ADAMS COAST GUARD STATION , 1955 
OF OBSERVATION AND SALINITY 0/007 
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On that days 26, lower low tide occurred at lower 


flow was 356,000 ots. 


‘Temperature 


as s indicated in in the tabulation: 
‘aye. TABLES 
_ TEMPERATURE OBSERVATIONS NEAR POINT ADAMS COAST GUARD STATION a 


14.8 


4 | 
22 
inity —— 
— 
 ap- 
ollege 
30" 
of temperatures at different depths were made by the Oregon 
_ State College and the Oregon Fish Commission off Point Adams Coast Guard 
r IF OSERUATION C TEMPERATURE (CEN) 


Januar 5 

ing the 1932 current survey indicated that the coarser bottom deposits were - 


found near the center of the section off Clatsop Spit and that seaward from 
~~ section the deposits became finer, the finest materials being found on the 


outer slope of the bar beyond the mouth. Bottom samples graded from coarse e 


_ to fine from the navigation channel toward the jetties inside and from fine to 


‘As the result of a study based on ey" of 199 —— taken in | the v vi- 


aa date. He believed that the materials derived from the erosion of the — 
ae lie on the ocean floor as a submerged coastal plain and that the beaches 


are constructional, built from sand that has come from the upper reaches of 
r 


a Columbia River and not from the adjacent cliffs. The clearness o 
freshness of the materials, their angularity, and the presence e of se several a 
particles adhering together were considered as proof of the transportation of 
the grains by floating during a recent period of time. He concluded that the 
_beach sands are not derived from the erosion of bedrock adjacent to the coast. 


but originate from a distant source. As the source was considered one of ial P 
about one-fourth andesite and two-thirds granite and rhyolite rock, it was be- € 
lieved 1 that the Cascade Gorge: could supply the former , but only the upper — a n 


reaches of the ‘Columbia River could supply the latter materials. u 


CLATSOP SPIT BIN SAMPLES 


DATE OF SAMPLING 


| 

7 
May 31, 1.55 
August 22, 1956 118 

August 22, 1956 2.00 
September 1957, “420°; | 


July 26, 
“August. 22, 


"September 16, 1957 


4 
August 22, 1956 
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of the opinion that the adjacent coast has in recent geological 
a 

- 

recent years s 
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_ Near the end of the 1956 dreding season an area of hard packed clay, ap- 

‘proximately 3000 feet long and 600 feet wide, was encountered at the 

level along the centerline of the channel near the turninalignment. ° 


‘ASCE COLUMBIA 
Referring to Fig. the alignment of the navigation down- 
stream from Harrington Point makes a long curve e to the south past Tongue 
7 Point and Astoria, thence through a a sharp reverse curve past Clatsop Spit | a | 


to enter the Pacific Ocean through the in a general south- 


cilities at Astoria and. vicinity, it appears ‘possible that a ‘change in alignment | 
so as to direct the major river flow due west from Harrington Point past _ 
Point Ellice. and thence to the entrance might be a factor which could bene- 
ficially influence « dredging - requirements « of the entrance channel and could 
also provide for greater ease in navigation. Such a modification in alignment 
=_— of course, require the maintenance of a channel connecting the re- 
~ aligned ‘main navigation channel with the | Astoria Harbor area and adjacent 


Preliminary of the Problem 
to Columbia. River, review of the aforementioned difficulties to navigation : 
points to the conclusion that there is a definite problem at this estuary __ 
_ entrance worthy « of detailed study in the interest of orderly and i timely plan- 
ning. In su such | problems, in which complexities are apparently ) compounded by 
_ unknown interrelationships of vaguely defined forces, preliminary analysis, _ 


based 1 upon knowledge and serves to lend sound di- 


- mation concerning conditions at the mouth of the Columbia River. Examination 
“of these surveys reveals that greater than project depths have prevailed since 
- the beginning of improvement throughout a natural channel, with length rang- 
ing | from 4 to’ 7 miles, lying just inside the bar between Clatsop Spit and Sand 7 
Island. Attainment of project depths across the aod in 1917 provided a \ direct 
connection between this deep channel and the sea, a connection which pre- 
-vailed until the formation of the inner bar between the jetties in 1952. Con- 
“currently, howeve er, the enlargement of Clatsop Spit pushed this channel - 
_ against the North Jetty and made it necessary to dredge the face of the spit tin 
_ order to maintain authorized project dimensions. As indicated previously, — i 
this: dredging, | confined to the Clatsop Spit area, totaled nearly 15, 000,000 = 
_ cubic yards during the period 1939-55. To obtain the increase in project depth 
to 48 feet, dredging since the a of 1998; has been necessary on both the < 
_ Comparisons of computed quantities of material - remaining within the 
authorized channel, shown in Table 1, with reported bin measurements of | 


_ quantities dredged indicate that recent dredging operations have | sometimes 


_ been accompanied by a natural scouring action which has greatly assisted the — 
_ Work of the dredges, and at other times by a heavy shoaling within the channel _ 


pears to be no consistent relationship of scour or shoaling to river adios, - 

| there does appear to be a tendency for t the > entrance to scour during a rising _ 
river condition and to shoal on a falling or almost stationary condition. Dur-— 
ing the 150-day period, between October 6, 1956 and March 4, 1957, when val ; 


river flow was quite uniform and about 140,000 cfs, ‘the 
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channel s shoaled 8.38 380, 0,000 cubic yards | as shown on Table 
This was partially compensated during the next 90-day period to June 3, 1957, 
when natural scour amounted to nearly 3,300 cubic yards. s. This shoaling 
trend was also evident during the winter of 1957-58 when the amount of ‘ma- ‘il : 

terial in the authorized channel was increased by a total of 5,850,000 + 
yards in the -day period between and February 19, 
_ The major portion of the materials dredged at the entrance saline 5 1956 . 

and 1957 has been dumped by the hopper dredges in an area lying about one 

mile south of the extreme end of the South Jetty. Condition surveys reveal B is, 
that, although some 14 million cubic yards of material have been Glaposed in 


that location since the spring of 1956, very | little of this material has re- 


Another phenomenon noted, particularly by condition surveys of April 1957, 


and later in November of that year, was a formation of waves of material, — 
ranging from 3 to 10 feet in 1 height, along the bottom of the entrance ve 
on both the inner and outer bars. The degree of wave formation within the i 
entrance channel area appeared nave pronounced than in any 0 of the adjacent - 
estuary bottom areas shown by these condition | surveys. 
ee As pointed out above, sieve analyses of bin samples of dredged materials 
_indic ate that these materials have the general fineness characteristic of sand, 
although on one occasion material coming within the classification of clay 
was encountered at the 48- -foot depth level. Previous analyses of the | compo- 7 
_ sition of these sands, revealing the presence of andesite, granite, and rhyo 
lite, led to an early assumption that their source lies in the upper reaches a 
the Columbia River. Although these upper reaches may have been the source 
. of a portion of the entrance sands during a previous geological period, the — 
| degree of shoaling occurring during the winters of 1956-57 and 1957-58 noted — : 
_above, when the upland discharge was relatively low, clearly indicates that < 
the upper Columbia River reaches do not now constitute the immediate source 7 
or locality fr from which material is now moved into the channel during, say, “< 
single year. _A possible source of these materials may be in the littoral drift — 
offshore along the coast; the evidence provided by the massive build-up of ma- 
_ terials behind the North "Jetty indicates that these littoral ‘movements, pre- 
dominately in a southward direction, may be of some significance particularly — 
when influenced by wave action. In addition to littoral drift, it is also possible _ 
that offshore bottom areas may be a source of the shoaling materials. An- 7 
other possibility, although seemingly less plausible, is that these materials x _ 
may have been deposited in the entrance area by wind action. Complete engi al 
neering consideration of the problem should, of course, definitely ascertain as Ff 


the immediate source of the shoaling materials. 
movement of materials within an estuary area, available knowledge nn 
long t been recognized as the primary agent responsible for movement of shoal - 
However, the action of river currents is known to be modified considerably ie 
_ A paper by H. B. Simmons(9) points out that engineers have been sguisigy y - 


these forces to river currents, density currents, littoral drift, wind, and i 
within an estuary by density currents created by differences in specific gravity : 
a 
in recent years that the amount of fresh into an sed) 


_ In considering the possible forces which may be responsible for the __ 
waves, or any possible combinations of these forces. . River currents have 

ing materials in nontidal streams and in the nontidal portions of tidal streams. 7 
of the fresh upland discharge and the salt water of the sea. 
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estuary and the degree to which it it mixes with the salt wa water of ee’ sea are — th 
major factors establishing the hydraulic and shoaling regimens of estuaries. 
With the source of fresh water at the upstream end of an estuary and the | a s 


_ i. source of heavier salt water at the seaward end, each type of water tends to | 
assume, as the result of density differences, a rough wedge shape with the — 
base of the wedge at the source, The sea water entering an estuary does so 
in this general wedge form, the entering tip of which moves along the bottom. 
> interface or line of demarcation between the salt and fresh water may _ 
‘from well defined to almost obscure depending on the degree of mixing 
of the salt and fresh waters. The ratio of mean upland discharge over a mean 
tidal: cycle interval t to the mean n tidal prism, as measured by t the volume of 
water flowing into an estuary from the sea during an average flood- tide pe peri- 
> od, may be used to roughly define three gee categories of mixing: highly 


in n the order of 0.25; and wit mixed, when the ratio is appreciably y less 


Distribution of current velocities in a highly stratified | estuary, such as i 
= Southwest Pass of the Mississippi River, shows that upstream from the 
tip of the salt-water wedge the direction of flow throughout the entire depth is 
toward the sea at all times with a \ a vertical | velocity ¢ distribution similar to that 
~ of an upland stream. Farther downstream the flow in the fresh-water strata + 
is still toward the sea, aes the direction of flow in the underlying salt water | 


ist upstream a at all times. The cause of the continuous upstream flow within 2 


q 


the salt- -water wedge is constant erosion of salt water from the interface by . 
the outflowing fresh water. The salt water thus eroded must be replaced to © 
maintain stability and, accordingly, a continuous upstream current is generat- 
ed beneath the ‘interface. In this type of estuary, rapid shoaling usually oc- 
curs in the region of the salt-water intrusion tip. The heavier particles of _ 
material, which are ‘Tolled or pushed along the river bed by the outgoing fresh 
Z water, come to rest as soon as the tip of intrusion is reached. ~The lighter 
_ particles, transported largely in suspension in the fresh water, gradually fall 


_ through the interface along the length « of intrusion and are transported up- 


stream by s salt water currents to the vicinity of the tip. -Flocculation of dis- ‘ 
_ solved materials in the fresh water may also occur along the interface and 
ce Available records indicate that upland discharge into the Columbia ia 
- estuary has ranged from an extreme low of 60,000 cfs to a maximum of about 
. 1,300,000 cfs. The ratios of these discharges to measurements of average 
= flood- tide flows shown in Table 2 indicate that the mixing of fresh and salt 
water within this estuary may assume characteristics of the highly- stratified 
estuary during tt the ‘spring freshets and of the partly- ‘mixed estuary during 
normal and low- -flow periods. _ Evidence that there exists pronounced salinity 
differences between the top and bottom layers of water of the estuary is avail- 
_ able from observations made by Oregon State . College and Oregon Fish Com- 4 
mission, shown in Tables 3 and 4, on July 25 and 26, 1955, at which time the 
upland discharge was about 360 000 cfs. This finding was confirmed by obser- 
vations of the Uz Ss. Hydrograph Office made in 1954 and 1958 as noted above. 
Earlier observations in 1936 indicated that the upstream penetration of the _ 
salt- water wedge may range from 11 to 15 miles. 
As the result of a recent preliminary reanalysis of the findings of 
~ Columbia River current survey of 1932, Simmons, in a later paper(10) stened 


_ that the effects of density are of major importance with respect to shoaling at 
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the Columbia River entrance. . While recognizing that more » detained iin 
later may show that certain of the resulting indications are no not completely 
valid, this preliminary reanalysis, made for three different upland discharges, — 
revealed that at a point 5 miles above the outer ends of the jetties, the flow 
at 90 per cent depth was predominantly upstream during an upland | discharge _ 7 
of 127,000 cfs and essentially in balance (50 per cent upstream and 50 per 
cent downstream) during upland discharges of 320,000 and 584,000 cfs. _Con-— 
sideration of only those current velocities in excess of | 0.6 feet per second, 7 
competent for movement of sand with median grain size of 0.4 to 0.5 mm (bin 
sample analyses noted above indicate median grain sizes of 0.25 mm on the — 
inner bar and 0.20 mm on the outer - bar) resulted in a marked increase in the 
upstream predominance of effective flow at 90 per cent depth for an upland 
discharge of 127,000 cfs, a change from a balanced flow condition to one of 7 - 
“upstream predominance for a discharge of 320,000 cfs, and no change from a_ 
balanced flow condition for a discharge of 584 000 cfs. _ Consideration of only — 
_ those current velocities in excess of 1.5 feet per second, at and above which = 7 
‘ane increase in 1 flood pesca at 90 per cent depth for upland dis- a q 
_ charges of 127,000 and 320,000 cfs, and a change from a balanced flow oe _—_ 
‘ dition to one of definite upstream predominance for a discharge of 584. ,000 cfs. 7 ; 
above available evidence and reanalysis of the 1932 current survey, to- 
_ gether with information not susceptible to complete check at this time, point . 
to great probability that density currents, created by salinity ¢ differences, 
constitute one of the major, if not the most dominant, forces controlling the ’ 
shoaling regimen of the Columbia River entrance. a: aciiiemaceimeaies! 
Except for its possible effect in creating the large build-up of materials 
behind the North Jetty, very little information is available concerning the — 
nature or extent of offshore littoral drift movements at the Columbia entrance. 
Considering the great mass of materials now trapped and consolidated behind © 
that structure, it may reasonably be assumed that an appreciable amount of _ 
material is moved by littoral forces across the entrance area. It may also 
_ be reasonable to expect that a portion of this material | may be carried a 
the entrance by flood tides and waves and deposited so as not to be — 
by ebb tides. Clatsop Spit and Sand Island might have had their origins in 
such action and study to determine the mineral and chemical composition | ona 
particle-shape characteristics, as well as mechanical analyses, should be 
made of the materials forming tl these geographic features to ascertain whether — 
or not they are formed of littoral materials. i 
7 fineness the existing project provides for the construction of spur Jetty 
” to reduce maintenance costs for dredging and to provide protection * 
against damage to the shorelines and main jetties from scour which may be ; 
caused by a change in currents resulting from the deeper channel, prelimi- _ 
nary analysis of the shoaling problem raises doubt concerning the effective- __ 


iD 


ness of that structure, particularly when viewed in the light of the probable _ 
forces controlling the shoaling regimen of the entrance area. Considering = 
spur jetty structure asa contraction work whose greatest effectiveness would <_ 4 
occur | during the high-flow periods, , it is doubted that Jetty “B” would be — 

instrumental in reducing, or affecting to any material degree, heavy shoaling _ 
Such as occurred during the relatively. low-flow w periods in the winters of eat 7 
—-1956- 57 and 1957-58. These factors, as well as as the cost of the structure, > 


) = estimated at $5, 944 000, and the possible economies which be 
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fore point to the the advisability of a more considered ap- Ww 

Preliminary analysis indicates ‘that such an n approach : should include a 


comprehensive s study of all phenomena susceptible of contributing to shoaling 
within the Columbia River entrance area, of which the phenomena of density 
currents and littoral drift appear to be the primary and controlling agents. 


Such a study should be directed toward answering the following questions: 


a, What and where are the immediate sources of shoaling materials? m 
= kh _ What are the effects of density currents and littoral drift on shoaling 7 : m™ 
c. If density currents and littoral drift are primary agents” for moving ma- al 

~ terials from the. sea into the estuary, what counteracting forces create a) 

the balance of conditions in the entrance pa ‘Pp 

In the light of experience gained in the y and of a 
entrances elsewhere in the United States, the problem at the Columbia River © i 
entrance was was presented ti to the ‘Committee on Tidal Hydraulics of the Corps of il 
Engineers early in 1956 and as the result of that presentation the Committee ‘a 
recommended that a program of prototype measurements be undertaken a ’ 
soon as possible with the view to investigating the problem in )a hydraulic — tl 
-modelatanearlydate. 


_ Available data which might be useful in a model study were examined, and 
_ based on this examination, , the | Director, Waterways Experiment Station, a 
Vicksburg, Mississippi in March 1957, ‘recommended a program of hydraulic 


and salinity measurements required for the desigu and verification of a hy- 


> 


obtain data for c of low, ‘normal, high river each. cycle 
to include observations of current direction and velocity for a continuous full 

; = tidal agra of about 25 hours at 19 stations located on six ranges as shown ~ 
- in Fig. 5. At each station measurements will be made on 30-minute intervals 


> 


at the surface and bottom and at the intervening quarter points of depth. 
Simultaneous observations will be made initially at one station on each rege: 
for a continuous period of about 25 hours to establish the relationship between 
ranges followed by simultaneous observations at all stations on each range 
until all observations have been obtained. . Salinities will be observed at hourly 
intervals concurrently with the observations of current velocities and di- 4 bs 
rections. practicable, an attempt will be made to obtain a few 
_ ments of current velocities, current directions, and salinities on a range be- 
_ tween the jetties. Considering the equipment which may be available for this” 
+ purpose, it is anticipated that a complete measurement cycle can be pti, 


it is planned to undertake the above program of prototype dindiaiiiiaaniti 
RB = as soon as necessary instrumentation, operating personnel and equipment, and 
can be made available. It is recognized that, even with an expedited 
(program of prototype measurements followed by the immediate 
struction, and verification of a hydraulic model, it may be a number of years 
4 before the forces that control the the regimen | of the Columbia River entrance iis 
may be accurately defined and suitable corrective measures designed and con- 


structed. In the meantime, it is anticipated that sufficient maintenance on 
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wil be made available to provide, by dredging, over 
a. 


An engineering approach to solution of the problem at the Columbia River 


“entrance should include a qualitative | definition of forces controlling the regi- 

_ men of the estuary area, such as may be obtained by means of a prototype | aoe 

" measurement program, together with a quantitative analysis of the magnitudes | 

and effects of these forces as may be provided by means sofa hydraulic model — 

; study. Such steps in the past have led to the solution of complex tidal 
problems in the Delaware River, Charleston Harbor, Savannah Harbor, and — 

_ other improved river estuaries by revealing the structural measures neces- 

sary to make effective use of the forces controlling the particular | estuarine 

regimens involved and similar steps should be no less instrumental in point-— 
ing to a solution of the Columbia River entrance problem. Even though as 
model investigation | might conclusively show that solution of this problem lies” 


tot in additional structural measures but in regular channel dredging, this 
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_ such a finding, substantial savings could be effected by Limiting future expendi- 7 
to maintenance dredging operations only. 
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- HY YDRAU ICS” DIVISION 
of the American: Society of Civil Enginece 
MODEL STUDIES OF HYDRAULIC DOWNPULL ON LARGE | GATES® 
P. Simmons, 


Air model studies were used to determine the hydraulic downpull forces — 

on a large fixed-wheel gate and on a slide gate. applicability of air tests” 
_ for this work is discussed, and the method of applying the data is detailed. — mm : 
= are presented that will be helpful in evaluating Gownpull rere aa oer ‘al 


The use of air models for studying | hydraulic pr problems has 
_ proved to be a reliable and accurate expedient. The advantages of air tests 7 
i ins place of hydraulic tests have been ' thoroughly discussed by Hunter —— 
and « others, and many examples of the use of these air models are given by 
J. W. Ball and D. W. ‘Appel. (2) Another unusual use of air models is pre- es 
sented in this paper, wherein studies o: of hydraulic downpull forces on a large — 
fixed-wheel gate and on a high- pressure slide gate are discussed. ee hy 
_ The downpull forces being considered here are those forces produced — a @ ‘ 
when 1 flow occurs beneath the gate leaf, thereby reducing the pressures | on +. 
the bottom of the leaf relative to the 5 pressures on the top. ‘This pressure A. - 
difference, acting on the cross-sectional area of the leaf, produces a down 
pull force which must be considered in addition to to frictional forces andthe _ 
leaf weight in designing the stem and hoist. The method used for obtaining Kt 
‘the downpull was to determine the unit pressures acting on the top and bot - 


areas t to compute the downpull fore forces. measurements of ‘stem loads \ were 


Note: Discussion open until June 1, 1959. To extend the closing date one ‘month, 
: written request must be filed with the Executive Secretary, ASCE. Paper 1903 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the 
- amertomn Society of Civil Engineers, Vol. 85, No. HY 1, January, 1959. yet 
a. Presented at the June 1958 ASCE Convention in Portland , Oregon. Bpetiat aa 


1. Hydr. Research Engr., ‘Hydr. Structures and Equipment Section, 
Lab., Bureau of Reclamation, Denver, 
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‘When the downpull first it was taken for granted 
that model studies would be made with water. . But it was ; later r recognized _ _ 
that air model studies would be equally satisfactory and that they would be © 
_ less expensive and time consuming. The method of f testing with air would be 
identical to that for water because the air velocities would be maintained be- - 
low about 250 feet per second, and incompressible flow equations could be = | 
used.(2) The degree of error would not exceed about 5 percent at a velocity — 
of 250 fps and would decrease rapidly with decreases in velocity. The _ 
velocities actually used were 64 fps or less, and the error was about 1 per-_ 
cent. Wetting would not be a problem in an air model and, therefore, un-— _ : 


treated wood and metal surfaces could be used. The construction could be 
relatively lightweight because the density of air is low and the weight and — 
pressure forces are small. The supply and exhaust systems posed no prob-— 
lems because the at mosphere served as both reservoir and receiver. _ 
eg argument in favor of hydraulic models was that by their use free — 

_ water surfaces could be. obtained. . With air models discharging into the at- - | 
= no free surfaces are possible. With either of the two gates under 
consideration free seamens will exist during free discharge operation. ae 


thought. In the case > of the fixed- marten gate, the 1 maximum . downpull occurs © | 
when there is no free surface and where back pressure or submerged condi- — 
tions exist. _ Thus, only tests with submerged flow were required. ' These are 
accurately reproduced in an air model. 
In the case of the slide gate, the downpull forces were desired at both free 
discharge and submerged conditions. After study of the leaf shape, it was | 
found that the pressure acting on the leaf could be analyzed for either free 7 
discharge or submerged conditions by proper treatment of the data obtained — . 
from an air model. el. The t treatment is discussed in detail later in this paper. 


UDO ws 


wale 
“The first was a fixed-wheel gate designed for the power con- 
"duit at Glendo Dam, Wyoming (Figure 1). This gate | has a leaf 3.3 feet thick, 
: 16.5 feet wide, 2 and 21.0 feet high. It is located near r the entrance of the 21- 
_ foot-diameter outlet and power tunnel and normally remains fully open. 
When a . regular closure is required, it it is 5 made w with no flow taking place me” 
with the pressures upstream and downstream from the leaf balanced. -Emer- 

gency closures with unbalanced pressures may also be made. Sustained 
operation w with unbalanced pressures will occur during the tunnel filling period 
when the gate will be used as the flow regulator at openings of 3 to 6 inches. - 
_ The control point on this fixed-wheel gate is at the upstream bottom edge - 
of the leaf. When the gate is is controlling the flow, the pressures on the leaf he 


bottom will be less than the pressures on the leaf top and less than in the tun a | 4 


nel just downstream. This occurs because as the flow passes beneath the leaf, 


the flow velocity— and hence the kinetic energy—becomes high. . Simultaneous- 


ly, and because the total energy in the water remains essentially constant, the 
The pressures on the leaf bottom were of greatest interest during the is 
nel filling period when the gate openings were very small. At these small © 
openings, the ratio of opening height to gate thickness was ‘small, so that le 


effect, leaf approximated a “short tube,” or more 1 
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specifically, a “short slot. ” As the flow passed beneath the — ‘it experi- 
7 enced it’s minimum cross section at the vena contracta, and then tended to 
bi expand and occupy most of the section between the conduit floor and the web — 
of the bottom beam at the downstream face of the leaf. ‘This action could pre- 
vent aeration between the leaf and the jet during free discharge operation, _ 
_ and would restrict the relief afforded by circulation of water during sub- 
- merged operation. Thus, lower pressures than usual would occur on the en- 
tire leaf bottom at very small gate _ openings, , and the downpull forces would ~ 
be extremely large. When the gate opening increases, the jet contraction on 
and subsequent expansion changes and more space is left between the jet and 
: the downstream web. The fluid is readily able to move upstream through the 
- _ space to relieve the low pressures on the gate bottom. When the gate reaches | 
the full open position, the pressures beneath the leaf are about the same as 


those in the tunnel and on top of the leaf. . At this point, the downpull is quite “i 
= “The second gate studied was the ‘Slide gate ‘design now being ea 
used by the Bureau of Reclamation. This gate is a rectangular, downstream _ 
seal, high capacity slide gate designed for use as a regulator under high Se et 
heads (Figure 2). The principal features of the gate are the small slots, the. 
a outwardly offset downstream slot corners followed by ‘gradually converging _ 
side walls, and a leaf with a flat upstream face and a 45° sloping bottom. This 
gate was developed to meet the requirements of the outlet works at Palisades - 
Dam where large flows of water at heads up to 240 feet are controlled. Db The Pi 
gate has since been included in the designs of several outlet structures with | *f 
heads ashighas 373 feet, 
magnitude of the downpull forces | on the fixed-wheel gate aS 
a on the slide gate were desired so that adequate, but not excessive, hoists _ 
and handling equipment could be provided. Model studies: were made on the 


Models 
air from m the atmosphere and forced it through a 10- -inch pipeline 6.4 sal 
of long, and then through the test section and back to the atmosphere a, 
‘(Figure 3). The rate of flow was measured by a flat plate orifice at the en- 


ons of the 12- inch- diameter r blower inlet line. All of the parts for "the ‘ya 


models were of simple a and inexpensive construction, and were fabricated in | 


The test section for the fixed- -wheel gate consisted of a plywood conduit 


16.67 inches high (Figure 4). The conduit was 8.25 inches wide upstream 
from the gate leaf, and 8.50 inches wide downstream from the leaf. _ The 2 a2- 


prooy on the bottom along or or near the conduit center line. Piezometers in n the 
tunnel roof upstream and downstream from the leaf were used to obtain the | 
pressure drop: across the gate. . The ‘slots were simplified and represented 
the width but not the | depth of the Glendo slots. The position of the gate leaf Ma 
4 could be observed through alae teaieane plastic windows that formed the outer 


i: The test section : for the slide gate was 16.67 inches high and 8.25 inches "7 a 
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Air was from the into inlet 
at center by centrifugal blower at left. 
Flow passed through a 10-inch-diameter pip 
section, then the test section, and back into 


3 - AIR AND 
and was aS supported in slots 0.64 inch wide. | ‘Three r rows of eight piezometers 
were placed on the leaf bottom so the pressure distribution could be deter- — 
mined at distances of 0.18 and 1.10 inches from m the sides of the leaf, and on on 
the center line. An additional ‘piezometer was ; placed on the bottom of one of, 
the projections of the leaf that extends into the gate slots. The upstream — 
pressure head was measured by a piezometer in the right wall 19 inches ae | 
stream from the leaf. A piezometer on the roof center line 1 inch upstream _ 
from the leaf gave the approximate bonnet pressure. No bonnet was included 
in the model, and there were no piezometers downstream from leaf. ie je 
_ All pressure measurements were made with water-filled U- tubes, and the { 
"readings were made in tenths of an inch and estimated to the nearest hun- ae 
_ dredth of an inch. The test procedure for each gate consisted of setting the 
leaf to the desired position, allowing a few minutes of operation for conditions 
“to stabilize, and then taking the pressure readings. ‘The leaf was then set to a 


nex ted. 
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4 measurement mr stem loads by strain gages, spring scales, weights, or othe 


devices. Previous experience at the Bureau of shown 


. Proper clearances between +a 

pa arts a are e of particular importance, and it is often difficult to predict pear “a 
Be. se clearances will be, or to produce them in n the model. _ When using the aa 


ee in » at regions ‘which ¢ can n effect ‘the downpull « on the leaf, and to a 


these pressures properly, 
In the case of the Glendo gate, the only condition where 


sure relief obtained beneath the leaf is that due to circulation of water be- a 
pressures on top of the leaf can vary from essentially 
_ when the backwater is insufficient to fill the tunnel, to high values when — 


‘aa no backwater is present and the gate opening is appreciable, ini sai | 

ill pass beneath the leaf and clear the downstream web enough to allow free 
aeration. ‘Thus, the pressures under the leaf will be about atmospheric. ‘The 

Brascnes oot above the leaf will also be about atmospheric because the gate 
‘seals are on the upstream face and the bottom surface is open to the » tunnel 
pressure at the downstream face (Figure 1). There is little or no downpull 

7 7 In the slide gate, downpull will be an important factor through a wide al 

_ of openings for both free discharge and submerged operation. The maximum 

_ downpull can occur with submerged operation because subatmospheric pres- 
sures can exist on the flat- bottomed portion of the leaf during this operation. a 
In contrast to this, the lowest pressures on the flat-bottomed portion of the _ 
leaf during free discharge operation are about atmospheric. The bonnet pres-— 


, sures will be | high at small gate openings for either free discharge or sub- 


_ merged operation because this type of gate seals on the downstream face of © 

the leaf, and is open to the tunnel pressures at the upstream face (Figure ~ 
The high bonnet pressures are the a factor in the high downpull forces 
produced on these gates. 


Wheel Gate 


Inthe simplified air ‘model values of head differential that will 
} occur across the full-size gate could not be set. But , as in a water model, 
‘the pressures on and across the leaf varied in accordance with the model 
eee and velocities. _ Thus, prototype pressures were obtained by multiplying 


across the gate to. the ‘measured model differential. To compute the protetype 
head drop, it was assumed that the 21-foot-diameter tunnel was ruptured at = 
the power bifurcation, and th that the reservoir was at the maximum water sur- 


face elevation of 4669.0, producing a head of 178 feet. The head loss — 

By losses become : smaller : and the head drop across the gate becomes larger. 

At very small openings the shut - off head of 178 feet is acting ~ 


—— 
: forces determined by stem load measurements because of the effects of fric- Bee ae 
— 
are cmall relative to the leaf thicknoece In this race the anlyv nresc- 7 
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January, 1959 
“the asi The relation of prototype head drop to gate opening is shown in | 


a 
The model was tested at a number of gate openings, and the pressure drop | Rv 


across the gate and the pressures on thi the bottom of the leaf were measured at 

each ¢ of the openings. These measured | pressures. are re shown in Table i 

interesting to note that the pressure just above the lower seal (Piezometer _ 

No. 1) is almost the same as the pressure acting on the roof at the leaf, an } 

"hence on the upper seal. ‘The v upper seal has the same plan area as the lowe 7 


one, and therefore the upward force on the upper seal is the same as the cell 


a 


downward force on the lower seal. These forces balance one another; thus =¥ 
-Piezometer No. 1 is ‘not used in the downpull calculation. = From these data, 7 


_ the prototype downpull forces were calculated as follows, using as an example ; ‘ 
Model pressures in inches of water we 7 


“Upstream Downstream _ web 
head head (average) ip 
average) 


Prototype head rop (from Figure 6A) = 178 feet 


“Under t the» web = 51) ) (21. 1.80) = = -11. 
the DS lip 52) (21.80) = 


 &Ifitis assumed that the pressure naan at any section along the leat 
is the same as the measured distribution at the center line section, the ‘down- 


pull may be found by multiplying the pressure difference between the gate top 
= bottom by the specific weight of water and ty the area of the = section 


_ > MODEL PRESSURES (IN INCHES OF WATER) ON GL ENDO F IXED-WHEEL G 


— 
— 
‘4 
es Downstream head = -(0.09 (21.80) = -1.96 = head ontopofleaf j 
‘Under the seal = -(1.25) (21.80) = -27.21 
Under the seal = -(1. .80) = -27. bi 
% 
Gate opening] 1 [ 2 [ 3 [ 4 5 [ 6 [| 7 | 
8.08| -0.72| -0.27| -0.26| -0.28|-0.27| -0.27|-0.26|-0.26|-0.26 | 8.06 | -0.04 
— 8.06) -1.25 | -0.50] -0.51}-0.53)-0.52] | 8.07 -0.09 
8.46] -0.97] -0.51| -0.48] -0.47]-0.49]-0.54] -0.50|-0.50]-0.51 | 8.53 | -0.07 
— 8.38] -0.81 | -0.53 -0.52] -0.56 | -0.56 -0.58 | -0.53)-0.50] -0.54] 8.58 -0.10 
7.57] -0.55 | -0.47] -0.46| -0.49]-0.48] -0.49| -0.48| -0.45] -0.49 | 7.74] -0.11 
7,14] -0.40 | -0.34] -0.34] -0.35 | -0.35] -0. 37] -0.36|-0.33] -0.41 | 7.28 -0.03 
-0.36 | -0.34] -0.35] -0.37]-0.37]-0.35| -0.34]-0.34| -0.35 | 9.20 | -0.05 
| 8.97] -0.05| 0.07] 0.06] 0.02] 0.04] 0.05) 0.08] 0.11/-0.01 | 9.36] 0.15— 
8.65| -0.66 -0.55| -0.55| -0.50| -0.47] -0.69 | 9.13 “0.49 
— 6.95} -1.92| -0.79| -1.80] -1.87] -1.85| -1. 86) -1.78| -1.72] -1.86 | 7.61 | -1.79 
— 5.56] 2.19] 2.28] 2.21] 2.23] 2.24 2.25| 2.25] 2.25] 2.25 5.60] 2.33 
66] 1.63] 1.67] 1.64] 1.64] 1.64 1.65] 1.62| 1.75] 1.63) | 


GATE OPENING-FEET 


= 


> 


- ES FEE 


RESERVOIR ELEVATION | 4669.0 


AGRO: 
FORGES ON ON GLENDO 165'x x 21.0! “FIXED-V WHEEL GATE 


DOWNSTREAM. 
HEAD DROP BREAK 1865 FEET DOWNS’ 
—_— 
4 | : 
DS. 
— 


1959 


96) (62.4) (16.5) (0.3802) 900 pounds 
1.96) (62.4) (08. 

“4. 96) 5) ) (0. 2240) = 


4 t 39, 670 pounds 


a A ‘aie te was made of the enman prototype downpull forces acting on the 
7 ‘seal, under the lower web, and under the downstream lip for a wide range of 

gate openings (Figure 6B). The total a. which is the sum of the ~~ naa 
> downpulls, is also shown. 
This downpull information r may y be applied other geometrically similar 


tallations b ll thod: 
installa lations by the following me metho 7a 


> Find the downpull that would occ ur at i red g e openings with the 
head drop for the proposed installation: 


= across gate, Figure 6A. 
3 Hy = Head drop for proposed installation 1 


B. pon the desired downpull by multiplying DP» by the ratios of the , 


areas under the two gates: 


under the proposed gate with | 
: under the tested gate 
= area under the gate 


_ bottom seal reached a dangerously low value of -27.21 feet of water. This is _ 
low enough to produce cavitation in installations at fairly high elevations. [fe 


Even at elevations as low as sea level, pressure » fluctuations of only moder- — 
_ ate proportions would be enough to lower the pressures momentarily to where 
__ eavitation would occur. If the 178 foot head differential across the gate were 


increased to 200 feet, the seal pressure ' ‘would become -30.60 feet and cavita- 
a could extend over much of the seal is At a 250-foot head differential, cavi- 
_ tation pressures would exist for gate openings from about 2 to about 4 inches. 


_ ‘For the conditions where cavitation is general, downpull on the seal is com-_ 
‘oo by applying ca aicaaiaes pressures to the area under the seal. Bo a 


2, 


= Downpull, Figure 6B 

DP 

Tt will eno ation tiat tie pressure beneath the 
iii 
i) 

— 


In the case of the slide gate, the flow patterns along the upstream > 
4 on > of the leaf and along the sloping bottom of the leaf will be the same _ 
_ free discharge and for submerged operation. Pressures measured on these 
components in an air model are therefore equally : applicable to free ee discharge 
and submerged conditions. A distinct flow difference will exist, . however, 
along the horizontal or flat portion of the leaf which lies downstream from the 
ny ending of the slope. During free > discharge th the flow will spring clear +e 
of this surface and the pressure on it, and on any seal extension, will be about 
_ atmospheric. During submerged operation an eddy of the flowing fluid will be 
in contact with the surface and the pressures will be a function of the head, 
the gate « opening, and the effective back pressure o1 on the gate. . The lowest 
_ pressure that can occur on the surface Ag water is flowing is the sll 
: pressure of water, about - -30 feet gate. ‘The pressure under a seal extension, i 


if one is used, will be affected in about the same manner SeTeeeY. i” au 


te - Practical considerations ry the width of the flat on the leaf bottom and 


may not be exactly similar for pice designed for different sized installations 
; and heads. To make the data applicable to these variations in design the data 

was reduced to that applicable to the sloping portion | of the leaf, and that ap- 
Tension to the flat-bottom portion, including the slot projections and seal ex- 


_—= bonnet pressures, and hence the pressures acting downward o on 1 top of 
, the leaf, were divided by the head differentials across the gate. This dimen- 
sionless ratio was plotted against percent gate opening in Figure 8B. ai 
larly, the pressures acting upward on the sloped portion of the leaf bottom 7 
4 = divided by the total heads producing flow. This ratio was also plotted in : 
Figure 8B. ‘The difference between these ratios, at any gate opening,  multi- 
plied by the prototype head the prototype differential 


this force, the forces acting on the flat portion of the leaf and on the slot 


projections and seal extensions are added. This sum equals the total down- ; 


The pressures acting on the sloping portion of the leaf were investigated in 
. Three rows or piezometers were included on the leaf, ‘one 0.18 inch a 
: from the leaf side, the second 1.10 inches from the leaf side, and the third in 
the leaf center line (Figure 5). The pressures for each row are plotted non-— ; 
_ dimensionally in Figure 7 for gate openings of 10 through 80 percent. The — A 
pressure at any point on a similarly shaped leaf bottom can be determined — 
_ by multiplying the appropriate dimensionless factor by the head producing - 
_ flow (Hy- hg) and subtracting this result from the total head upstream from — 
"the gate. For determining the downpull forces, the model pressures were — 
~ replotted dimensionally with the forces being considered vertical. Planimeter . 
measurements were made of the areas within the pressure distribution en- 
vce = for the sloping portions of the leaf. These areas were divided by the 


cavita-_ 
al, cavi- 
inches. | 


- > .23 inch thickness of the leaf above the sloped bottom to give an average 


pressure value at each piezometer row. 4 The average pressures obtained for r 4 
_ the row 0.18 inch from the leaf side were assumed to act from the leaf sides _ F 
2 stations 0.50 inch from the sides. The pressures for the row 1.10 inches ~ a 

from the side were assumed to act from the 0.50 inch stations to 


= 


y, 1959 
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1 
_ leaf. This head difference, multiplied by the appropriate area, and by the a | = an 
< ensi a ate : in the downn thi na rt he le = x 
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HYDRAULIC DOWNPULL 


“inches from the sides. The pressures at the center line siaieibiaas were ce 
assumed to act over the remaining width of the leaf. The products of these © 
“average pressures and appropriate areas were added together for each a 


- opening. These pressure-area summations were divided by the projected 
"area of the sloping portion of the leaf to give the equivalent single pressure 
a values 3s that c can be assumed to act over the sloped area of the leaf at each © 

"gate opening. These | pressure * values, divided by the total head producing — 
— flow , establish the curve plotted in Figure 8B. At zero opening, the average r 

_ pressure acting upward on the plan area was the same as the total head. It z 
decreased to a minimum at 70 percent open and increased again at 80 percent 
open. No tests were at above 80 percent because the d downpull 


would be omen. 


= 


In the full- sized slide ines ain tieen of the leaf carries a seal that is in | 

contact at all gate | openings: with a seal plate in the bonnet. . There can be no o 
4 flow or appreciable leakage between the back of the leaf and the bonnet, and 

only a little flow down the gate slots. _ Thus the bonnet pressure is about the 
= as the stagnation pr pressure just ahead of the leaf. The model pressure 


ahead of the leaf showed that without submergence the bonnet pressure was _ 


the same as the total head in the conduit at zero gate opening, and that it de- 


: measurements, based on the piezometer in the top of the conduit linch 


__ creased with respect to the total head as the gate was opened (Figure 8B). | 


‘The difference between the bonnet pressure and the average pressure over 
the sloped bottom was a maximum at an opening of about 45 percent and the 
difference did not exceed about 53 percent of the total head (Figures 8Band © 
C). The net downpull on the flat bottom, , the slot projections, and any seal 
extensions can be determined from their areas and the difference between 
the bonnet pressures and atmospheric pressure. The total downpull will be * 
the sum of the individual downpull forces acting on the various asad of the 
: A sample calculation of the downpull forces follows: 


Assume: discharge conditions, gate leaf 6.0 feet wide, 2.0 thick, 
3- -inch flat on bottom, 6- by 6-inch slot projections, no seal 
tension (Figure 8D) Head H = 200 feet. 


acting downward on the plan area over the sloped part of the 
im” 23.9| 4 52.9 


42.6 94.2] 105.0 105.8] 102. 


downpull, Ib 
(Area = 6.0 x 
1.7! = 10.50) 


q 
> 
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| 
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Bonnet pressure 


Maximum downpull on this’. “Downpul on this portion is 
portion is 53% of head free_dischorge-bonnet 
producing flow times ‘head minus atmospheric 
pian area, times 62. head, times area, times 624 a 4 
A METHOD OF COMPUTING Submerged discharge - -bonnet B. PRESSURE IN BONNET AND ON 


head minus static head 


seat, times orea, times 624 


GATE OPENING-% OF EFFECTIVE TRAVEL 


Sate width 60 
Gate thickness 2 0° 
Flat bottom width 3° 


Gare OPENING - -% OF EFFECTIVE TRAVEL ATE OF EFFECTIVE TRAVEL 


NET DOWNPULL PRESSURE FOR 


AREA 


& 
7 
ESSURES IN BONNET AND ON LEAF, 
FIGURE CAL DOWNPULL CURVE FOR 
— 


_ HYDRAULIC DOWNPULL 

cal The: net force acting downward on the plan area of the flat bottom and the - 

“glot projections (no seal extension on standard gate) is: ghee 


Bonnet Pres 


Bonnet Pres im ve" 


Head diff (ft, 
HAO) 
62. 
(8)x Area= [24 
Downpull; A= 
(0.25 x 6. i 
2(0.5 x 0.5) = 


aw it were assumed that the full total head acted over the entire sectional 


= = 156,000 pounds. . The maximum downpull of 91 ‘810 pounds obtained with the 
. _ above test information is 59 percent of this value. This percentage will vary | 
- somewhat depending upon the width of f the flat bottom of the leaf, the location 
of the leaf projections, and whether or or not seal extensions are used. aes nea a 
oe If the tailwater elevation rises to submerge the gate or fill the downstream og 
conduit, and the headwater elevation remains unchanged, the head differential 
across the gate, and hence the rate of flow, decreases. In spite of these | 
a. and ii that the Reynold’ s Number remains reasonably high, 


: ‘difference between the upstream conduit total head and the head at any ; point ‘di 4 


on these areas of the leaf, divided by the head differential across the Pes | 4 


remains the same for any gate opening. 


K, and the downpull force on the sloped portion of the leaf will be q 
the same as for the free discharge condition with the same differential head. — 
‘ the case where a gate can be placed more deeply or less deeply below the _ 
-tailwater elevation, the effect of submergence change is merely to equally 


raise or lower the bonnet and leaf slope pressures. _The head mu orces 


no cavitation at the lower submergences). 


When the discharge is into a pool, the head producing flow may be taken as : 


> difference between the total head upstream from the gate and the pool — _ 


Z depth downstream. When the discharge is into a filled conduit, the head 


2 producing flow may be taken as the difference between the total head 


— 
| 
| - 
| 
| 
the leaf ramaing the eama And heraues the flaw the game the 
| 
| 
q 
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- “January, 1959 AS 
For submerged flow where conditions would permit vapor pressure under 

the flat on the leaf, the method of calculation would be the same as that in the 7 
a free discharge example except that the downpull head differential in Item 7 _ i. 
-would become the total bonnet pressure plus the numerical value of the vapor : ; 

_ pressure, about 30 feet of water. This produces the greatest downpull pos- 

sible, | but it is unlikely that vapor pressure would ever occur over this whole 2. 
Thus, this extreme downpull condition will probably never be reached. 
In the more likely cases where conduit failures or unexpectedly low tailwater i 


will not t be fac tors, , and where the eeeageaee will always be sufficient to 7 


net downpull would be determined from the difference between the bonnet 
‘ pressure and these surface pressures. Such surface pressures may be diffi- 
— cult to predict exactly, | and would probably have to be determined experimen- 
tally if precise results were needed. But for ordinary work, the 
can be obtained from the non-dimensional data in 


SUMMARY 


Air models are reliable tools for flow phenomena and pres-— 


submerged. _ They may also be used for downpull studies when the gates dis- a 
ree freely, provided that care is taken in analyzing the flow conditions and 
‘resulting | pressure distributions on the various portions of the gate leaves. * 
Air models are less expensive to construct and test than hydraulic models, 
and they are accurate, convenient, and easy to work with 
_ The two gate types discussed in this paper are typical of installations being 
“made throughout the country. data are presemes in a form usable for 
_ geometrically similar to the test gates. The maximum . downpull c on the ray 
- fixed- wheel gate occurred during submerged operation at about a 3-inch open - 
ing. The downpull rapidly decreased as the opening increased. On the slide 
gate, the maximum downpull force occurred at about a 45. percent gate open- o 
ing. This calculation assumed that the total head remained constant at the © 
gate regardless of gate opening. In the usual case, friction in the ‘conduit: 
: ahead of the gate w will cause a ‘loss of head as the flow rate increases. . The l» 
maximum downpull would be less than in the example, and would occur at a y 
"smaller gate opening. i Bonnet pressures play a particularly important part in 
the: case of the slide gate, , and in designs where leakage from the bonnet is 


ee The studies discussed in this paper were conducted in the Hydraulic 
Laboratory of the Bureau of Reclamation, and the fixed- at gate studies 
_ were made by Donald Colgate, Hydraulic Engineer. pas? n 
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> vapor = prem oi Bulletin No. 31, Studies of Engineering, State University of lowa, a 
pos- City, 1947 (Proceedings, Third Hydraulic Conference). 
Pray 2. . Ball, J. Ww. , “Model Tests Using Low Velocity Air, with discussion by 


Paper No. 2517, Transactions, Volume 117, 1952. 
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APPLICATION OF SNOW HYDROLOGY TO: THE COLUMBIA BASIN» = 
liver -Johnson,! M. ASCE and Peter B. 


Observations of the weather elements: cause the of a snow 
pack were made for several years in three laboratory areas located in dif i 


i ferent weather regimes of the high mountains of the western United States, ~ 


and the data obtained were analyzed | and correlated. Starting with the inter-_ 


| relationships between the weather elements that produce snowmelt at a point, 
| methods and procedures were found for the rational analysis of snowmelt in “4 e 


‘relatively small basin. The application will be enlarged to , cover large parts on 
| of the Columbia Basin as the network of stations covering the basin is en- | 
larged ar phenomena are measured. 


in the belt of the o heavily woistere- laden prevailing westerly winds. In winter | 


the Aleutian lows attract maritime air masses which deposit their moisture _ 
as they ascend the high mountain slopes on the west side of the Cascades. The 


in n the “ rain shadow” of the Cascades, and annual precipitation in n this: region * 


may be as little as eight inches. Additional uplift of the air masses by in- mod 
terior mountain ranges ai and the Rocky Mountains causes increased oe 


tion in the headwater regions, in amounts ranging to more than 80 inches per 


Note: Discuss Discussion open until June 1, 1959. To extend the closing date one month, a 
a request must be filed with the Executive Secretary, ASCE. Paper 1904 is 
part of the copyrighted Journal of the Hydraulics Division, Proceedings of the viet = 
of Civil ‘Engineers, Vol. 85, No. HY 1, January, 1959. 


1. Hydr. Engr., “North Pacific Div. ., Corps of Engrs., ‘Dept. of the Army, 


2. Hydr. Engr., Portland Dist., Corps of of Engrs., ‘U.S. Dey Dept. . of the Army, 


Se 
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d 
Columbia Rix geannual 
runoff ranks very high among the streams Of the United states, averaging = 
id 
— 
| 


‘The Columbia rises in Lake at elevation 2652 feet m.s.l. 


- Compared | to elevation 1670 for Lake Itasca, the source of the ee 


= 7 The mean elevation of the Columbia . drainage east of the Cascades is Canal 


_ 4500 feet m.s.l. The rugged topography is important because of the large > hi 
-vertical al distance through whi which the water descends on its way back to the | Sl 
; Pacific, and because of the possibilities for yr economical construction of large | W 
dams afforded thereby. In January 1957, 25 million kilowatts of hydro capaci- on 
_ ty were available in the United States while that in Oregon, Washington, Idaho CE 
and that part of Montana west of the Continental Divide (comprising about 9 -.. 

_ percent of the U. S. area) totalled more than 7 million kilowatts or 28 percent § t 
— During the ‘ites of a storage project on the Columbia or one of its tribu- # 

_ -taries , snow hydrology must be considered in the computations required in the dl 
selection of the project design flood. This is the largest flood that must be 
controlled if a satisfactory flood control plan is to be achieved. _ For that part § ° 

_ of the Columbia and its tributaries lying to the east of the Cascade Range, the ( 
melting of the snow pack creates the severest ‘design conditions. 3. In the Willa- 
mette and other coastal streams the most intense floods are caused by none 7 
rain- on-snow events occurring in the period November through February. | P ( 

és _ Power projects and flood control storage are entirely compatible in the s 

; _ storage projects above Bonneville Dam. Power draft through the fall and i 


early winter will provide ample space for storage of floodwaters from moder- 
ate rainstorms that have occurred infrequently | in December, and saan. 


_— draft w will provide ; additional space before the onset of the snowmelt season. 
oe _ In the storage projects above Bonneville Dam there are monetary a 
to be derived by releasing stored ¥ water for power production at site and at 
downstream | plants, and there are monetary benefits derived from flood -_ i 
damage prevented by such storage, both locally and in the vulnerable reach: 
of | the Columbia lying downstream from Bonneville Dam. — 
a In the operation of the Columbia River | projects, the p principal eis is 
to provide storage space for a possible snowmelt flood and still have full 
reservoirs at the beginning of the power draft period. i To achieve this ob- f 
jective, the water content of the snow on the ground is measured at frequent 


intervals throughout the accumulation and depletion period. With the current | 
: _ observations of precipitation and temperature and assuming normal weather 
for the remaining period, satisfactory runoff volume forecasts are made. - 
a These forecasts are increasingly more accurate as the season progresses © . 
_ because t. the future precipitation becomes progressively | less important. Large 
"deviations from normal precipitation during the melt period in May and June 
can cause the runoff volume to differ considerably from the forecasted event, 
; but operating parameters can be devised that w will afford satisfactory r sons 
_ There are other snow-fed streams that are of great importance to the a 
United States, such as the Colorado, Sacramento, and Rio Grande, to name a 
a few on which the Corps of Engineers and Bureau of Reclamation have built | 
many reservoirs. JA knowledge of the rates of accumulation and depletion of 


‘as also have a great interest because they are the principal agencies 
_ charged with th the collection, processing and dissemination of hydrologic in- = 
7 formation. All are endeavoring to improve their forecast of the available 4 


basin snow water, of its rate of release and of losses. / Reliable forecasts by 
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SNOW HYDROLOGY 
- relatively | rapid means is mandatory for both economical design and o oper 
See need for detailed, accurate knowledge of the characteristic and be- 
_ havior of snow was ‘recognized by several of the Federal agencies and after 
- some discussion, a program m known as the “Cooperative Snow Investigations” 
was agreed upon and supported by the Weather Bureau and the Corps of Engi- 
_neers in 1945. Limited assistance was furnished by the Geological — 
H Bureau of Reclamation and Forest Service. Direct participation by the “eg 
i Weather Bureau ended in 1952, and thereafter the basic research was con- 
tinued by t the Corps of Engineers under the name of “Snow Investigations.” ina : 
‘Throughout its life the program was pointed toward the solution of a 
problems pertinent to mountain regions of western United States. “Three 
laboratory areas were heavily instrumented and studied intensively for peri- 4 
ods of five to eight years in order to provide the basic data for determination 
a the hydrologic factors involved. These laboratory areas were designated - 
Central Sierra (near Soda Springs, California), Upper Columbia (near Marias 
Pass, Montana), and Willamette Basin (on Blue River, in . Oregon). enact 
| ieielae in 1954 a part of the group effort on snow investigations was : 
‘directed to toward d compiling the results achieved | by these studies in the form — ‘od 
of a summary report that would serve as a reference volume for hydrologists. _ 
The work was concluded in June 1956 with a volume titled “Snow te og: il 


: Ber. The field c of snow hydrology is concerned: with the evaluation of various - 
components: of the hydrologic cycle. It draws upon the sciences of meteorolo- . 


and Areal 


gy, hydraulics, thermodynamics, geology, soil mechanics and plant physiology. 5 
‘The development of the applied snow hydrology begins with observations and _ 
-measurements, at a point, of factors affecting each component. This is fol- 
lowed by an analysis, leading to the establishment of mathematical relation- _ 
ships. Prior to basin wide application of the point values, the relative magni- 
tudes of the components are determined, and emphasis is placed on the more an 
important elements. Having determined point values of the components of the 7 

hydrologic cycle, the next problem is to modify the equations to fit the local — = > 
climate and topography peculiar | to an individual project basin. During the 

af course of the snow investigations it was discovered that experimental data r 
from areas as large as the laboratories (4 to 20 sq. mi.) could be fitted to | . 
= relationships only with difficulty, and near the end of the program =< 


perimental data from impervious areas measured in the hundreds of square © 
feet were obtained to assist in evaluating the coefficients needed to determine 
: _ Having deter mined point values, the next step consists of utilizing these va 
values to determine amounts and distribution over a basin area, with respect 
to differences in environment, as well as to time. This involves procedures : 
which are much less exact than those for point evaluations. _ It is impractical | 
_ to attempt a point-by-point analysis in basin application; rather, it is neces- 
sary to d deal with basin averages in major subdivisions of geography or en- 
vironment, and also to deal with averages in time. This concept leads to the 
se of point measurements as indexes to represent basin 


— 
the 
d 

ade. 
— 
ea 
built 
tion of &g 
| 
A 
to 4 


a> 
> 
use _— indexes involves the theory of sampling errors in ‘measurement, statis- 
_ tical procedures may be employed to evaluate the reliability of estimates and — 
also | the e weights assigned to individual factors. ” Intelligent selection of indexes 


requires an —— of the heat processes and the relative effects of 4 
ard 


‘Water - in Areas of Snow Accumulation 
reliable are available, but the records too short 
to to justify a statistical approach, an accounting method, called water- balance, = 
_ is employed in forecasting the stream runoff f during a Specified period. 4 = 
addition, a knowledge of the water-balance for a given area is necessary in 
order to select appropriate forecasting parameters and to interpret their re- 
liability. The water -balance technique may be used as a forecasting proce- ' 
_ dure by quantatively balancing runoff against precipitation, ‘change in snow- 
pack water equivalent, and losses. A knowledge of the water balance is also 
necessary for rate-of- -flow forecasting and for design flood tae ened 


The water- balance method is formulated by the equation: 


in which P is on net ong iain (after correcting for gage catch inn 
ficiency and interception loss); W, and W2 are the basin snowpack water — 


equivalents at the | beginning and end of the period; and L is the loss « con- 


sisting of evapotranspiration Maa) ant soil moisture change (Ls), that is, 


in | which Bis is positive when water is a added to the soil to redace its moisture 
deficiency < and negative when | moisture is being evaporated ed from t the soil; Qe 

is the generated runoff volume during the period considered, equal to the aa 
served runoff corrected for recession flow. Generally the purpose of a run- 
off forecast is | to determine | in advance the volume of water (Q) expected | to 3 


B: which Q; and Q2 are the runoff volumes under the initial and terminal 4 


oe is isbn to ‘melt the snow. Atp present we know ‘enough about the jae d 
; principles of heat exchange between the snow and its environment to formulate 7 
snowmelt dh based upon detailed studies for each of the three 
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» SUE the heat processes. The level of our confidence is ; not as high in in the enatines 
fes and — tion of the point snowmelt principles to basin-wide melt, ‘since only a small 

indexes umber of basin-wide melt analyses have been completed. basin pre- 
Ss of cents a separate snowmelt problem and has a unique solution. The nsec 
B esig is to integrate the effect of topography on the climate near the snow. 
This involves (1) the separation of the project basin into areas of homogeneous 
-} heat supply under cloudless, cloudy, and rainy sky conditions, and (2) deter- — 
“mination of the basin constants for each area—constants which will integrate 2 
- the influence of topography on solar radiation, on long wave radiation, and on 


Pri fis Ther processes F responsible for the supply or less of heat are: > 

eir a | Net longwave radiation exchange between the snowpack and its environ- 


Convective heat transfer (sensible heat) from the air (He) 
Latent heat of released by condensate (He) 


Conduction of heat from underlying ground (Hg) 
content of rain water (Hp) 


“These processes are influenced by 1 ‘meteorologic, ‘geographic and time factors 


which must be considered and related to the ordinarily observed data. The es. 
= equation for snowmelt at a point is: igh 


which M is in inches, is the sum of heat energy short 
moisture | > wave radiation, long wave radiation, convection, condensation, conduction and 
oil; Qg | Precipitation; 203 (= 80 x 2.54) is the number of calories required to melt one 


the ob-— inch of water equivalent of ice at 0° C over a unit area; Bisthe thermal | 
arun- | quality which is equal to the ratio of heat required to melt a unit weight of the 
ted ito snow to that of ice. B>1 when the snow temperature is less than 0°C and 


erefore 1 when the snow is at 0°C and dry. Buta melting snow, after 

- contains liquid water which requires no heat to melt. _ Therefore, its thermal © 
ae "quality, after drainage, varies somewhat with density. An average value me. 

iM = 0.97 is satisfactory, as shown in equation (5). ae 

ia The analysis of observations and measurements of radiation, temperature, 

‘ humidity, condensation (or evaporation) and melt at a point resulted in the 


following equations, where snowmelt, M, for each heat a is in inches | 
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is the incident shortwave radiation in langleys (Cal. /sq.cm.). 


Longwave Ra ye Radiation we 


‘For clear sh skies = 

is the air temperature at the 10- foot level in OF, 


q where 


+ For cloudy (low 
M 0. 029 (T, 32) 

eX, is the he temperature of the cloud base ‘in oF 


Convection BM Melt (th 7 


is s the ratio of the atmospheric pressure at the site to ans 


where p/p, Po 


Za is the height in feet above the snow surface of the air temperature 


the neigh is in feet eet above the snow surface “ the 1 wind ‘speed measure- 


the air in OF. 
_ is the snow surface | temperature. 


i. is the wind speed in miles per hour. 


Condensation Melt It (theoretical) 


maid 


where’ *. is me height in feet | above the the snow ‘surface of the wap vapor pressure 


Zp and are as defined above. 


1 = is the vapor pressure of the air in millibars. 
eg is the saturated vapor pressure of the snow surface. — 
Combined Convection and Condensation Melt (simplified) a 
= (70.008 v7 22(Ta - 32) 0. 78 (Tq - 


where v is is the’ e’wind speed at the 50-foot level in ‘miles p per hour. 
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ASCE 


of Heat from Ground 

_ Melt from this source is negligbly emelt, so that direct evaluation is un- 


“necessary. Average value is ‘approximately 0.02 inches/day. 
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practical method of a reliable clear -weather melt equation, 


E in terms of the available heat indexes, for a ripe snow in a given basin is by 

- means of statistics. The influence of topography and omen cover on each © 
= process will then be integrated in derived coefficients that will also in-— 7 
Fi clude the period of melt and the constants of proportionality between units of © 
measure. proposed form of such a linear regression equationis: 


where is the generated runoff. The term, k! (1-a)lIj, is the solar radia- 
_tion melt in which ]j is the solar radiation incident on the horizontal surface ; 
of the snow area, a is the albedo of the snow, and k' is the basin constant 
integrating the effect of topography and forest cover on solar radiation. The — 
term CT, defines the long wave radiation melt with C as its derived coef 


ficient. term, Kee(BT a+ CTq - Tp) v, is the convection -condensation 
melt, in which = is a constant defining the degree of basin exposure to 


wind; T, and Tg are the air and dew point temperatures; B, C, and Tb are 
derived constants, and v is the wind speed. D is the regression constant. ce 


= melt is negligibly small and it will be partly ental in the other 


The melt at the snow surface one the generation poring Se: is: rede 


in which L is the loss of water by evapotranspiration and deep percolation. 

Note that the proposed statistical analysis assumes L as a function of > 2a 
and I. Evapotranspiration is a function of these elements; and deep percola- 
tion, unless definitely established, is assumed to be negligible. il eee, 

ao Examination of data taken on « cloudy and cloudless days indicates that an 
_ adjustment factor for average cloudiness should be used in computing long 
wave radiation melt. On the ‘days of complete cloud cover and on rainy days, 
it is easier to develop surface melt equations. by a rational approach, 
j “ate equations peculiar to each sub-basin may be derived statistically by 
direct correlation of appropriate meteorologic f factors to daily generated run- $ 

“off amounts as s explained above. In many cases, it may not be feasible to de- 
: rive basin snowmelt equations in this manner because of lack of appropriate _ 
data. General snowmelt equations were therefore developed from laboratory 


% data, which may be used to evaluate basin snowmelt on the basis of the local 


— 

where T,, is the in (usually assumed to be that 
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— 

neasure- 
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— 

physical characteristics of a basin affecting heat supply may be expressed in 


> 


.* solar radiation. The form of the equation is determined by the amount of © 
ae ‘forest cover, since it has a marked effect upon the relative importance of ae 
both wind and solar radiation in estimating snowmelt. . For rain-free — _ 


y > modified by the appropriate basin factor to account for slope and orientation — 
(k'), forest cover shading (F), to wind id (k), and for cloud cover 


“resulting in the following equations: 


Shortwave radiation melt, Mrs =" —* whe 


— | 


N) (1-F) .022 -32) 32) - 
= 
0.0084 2 (7 0. 78(Tq -32) 7 


In the above equations Ij is the daily short wave radiation (in sali " 

‘reaching 2 a horizontal surface at the latitude and altitude of the snow area; a 
Ta and Tg are the daily average air temperature and humidity at 10 feet — 

above the snow; ana V59 is the daily average wind speed (in mph) 50 feet —_— 
above the snow. This study was followed by a statistical analysis of the clear 
_ weather snowmelt in each of the three laboratories, using estimated basin 4 

average temperature, wind, k, and F values. The analysis | resulted inthe 

following typical | snowmelt formulas, in 1 which k, , k! an and F are estimated from 


a careful examination of the topography of the snow EEE ER " 


forested area a (like Willamette Basin ‘Snow Laboratory 


M 


=k (0 0.0084 v v) (0. 227, 0.78 Tq) 


wind speed, solar radiation and cloud cover. During periods of active melt, 
only variable of the snowpack condition is the albedo and the snow 
For estimating snowmelt during rain, values of air temperature, wind speed, 
ae 
| 
— 
at > 
4 
i 
ry) having 25 to 
EEE 


1959 | 3} 
(0.0051I, j) (lea) -N) (0. .022 Ty - 0. 84) 0297, 


(0.0084 v) (0.227) # 0.787; >. 


+0084 
wheres 


MO 4 is the snowmelt in inches per day over the snow covered area 


is the difference in air temperature measured at 10 feet above e the 


‘snow and the melting snow surface temperature, in’ OF, The fore- 
ms - casted temperature of the air at a station may be extrapolated to the > 
snow level by the lapse rate. The of the is 


is the difference b between the -dewpoint temperature at 10 feet above 
La the snow and the n melting snow surface, in OF. —iItis approximated 
_ from upper air or surface humidity forecast. _ The dewpoint of the 


| melting snow is 32°F 


is the difference between the cloud base and the. melting ‘snow surface - 


7 S temperatures in OF. It is estimated from the upper air temperature - 
or. by a lapse rate from the surface temperature. 
vO is an estimated a -ssprpuecacwsnhiad wind speed at 50 feet above the snow, in 
miles per hour. Wind stations in the proximity may serve as a 


i, 3 the estimated short wave radiation on a horizontal surface at the the 
elevation and latitude of the snow area in n langleys perd day. It may be 


from duration of sunshine, cloud cover, or diurnal temperature varia 


a _ 4 mated from more 1. A melting snow has a smaller albedo than the “a — 


aes and falls very rapidly to lower values with time, following a decay- ey 
el _ type curve. _ The accumulation of residual impurities, as snow melts, 
tends decrease the albedo and in it be as low as 0. 


is an ‘estimated average e basin cover, effective in in shading the area 
solar ‘radiation, estimated asa decimal fraction. The estimate 
5 im of F involves the consideration of the type, — and density of the | 


is a factor intended to correct Ij for an average slope angle and 
_-—-—s orientation of the snow covered area. The effective short wave radia- zi 
dec decreases as a sloped surface is oriented from south toward 


north, and as the slope angle is farther away from normal to the sun’ s 
rays. effect k' = I/lj, the ratio of the insolation on the true snow 
_ ‘surface | to that on its horizontal projection. . Figure 2 shows I/]j values — 


fort the entire year for a southerly and for a northerly oriented unit a 
area of 25° * slopes in 4¢ 46° 30 N Latitude at 5200 feet above mean sea 


ibe 


ko ot an average corrective factor for the | degree of exposure of the snow a 
area to wind. This factor has a ee at 1. 0, , but it atid — 
as low as 0. 3 for a forested = 
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is the estimated cloud cover, _ expressed as a as a decimal fraction. _ 


to equa » Vy and k are missing. 
a forested area, the effect of wind is to a near constant 


melt runoff, essentially includes convection-condensation and long \ wave os 
tion. These components are satisfactorily defined by a linear function of the © 
air temperature (T) and the dewpoint (Tg), within the range experienced poe 


proach resulted in the coefficients shown in i equation | (19) for for the — 


ing the melt period. AC combination of both the rational and statistical a? 3 


proach , assuming insolation melt plus ‘ground melt equal to evapotranspiration; 
and therefore, the generated daily snowmelt runoff (M,) equals the convection- 

__ When equations (19) and (20) are used, “water excess” for routing is ob- | 

_ tained by subtracting only the infiltration water which is not available to direct 7 

runoff within the hydrograph period. Equations (21) and (22) are derived by | 
rational approach, considering point sarface melt. The snowmelt M is the 
water input at the surface and includes the evapotranspiration. 

_ The snowmelt during rain was derived from clear weather melt equations. 

Since e Tq = = Ta isa permissible ‘Substitution during a rain~ -on-snow event, the 


daily c convection- condensation melt from equat mn (18) for an open to part 
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RADIATION ON A “HORIZONTAL SURFACE 


Assuming the temperature of the cloud base, forest foliage, and the air 


the rain melt a 


_ The short wave melt varies with the time of the rain. But it is very nae 2 


Small and it is considered satisfactory to to engress | it as a constant or as a ee 


———| Direct solor radiation only 
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areas. 


10 ion lightly forested |and areas, or 


melt is considered if any, a time: when the 

heat supply may be used up by the moving surface and subsurface water. Melt 

The components are combined, and the following « equations are written to 

: represent the daily snowmelt values during a rain-on- snow event: es } 

to forested areas, (23), (25), (26) a and ind (27) 


M' - 32) (0. 0. 007 oO. 0. 05 (2-F F) (28) 


For forested to heavily forested areas, using equations (24), (25), (26) and 
. 
Note for a k and F are constants. "When they are e determined, 
equations are further simplified, allowing the construction of a 
table or nomograph for zonal or basin melt values. 
the entire basin | (Mp) is: 


wh re is th of the snow covered area | to the area of the entire 


a Effect of Snow Pack Conditions on Runoff 


: 


it may be retarded 0 or detained by - ice e dams, by ‘choking, and by suspension ‘oi 
e a form of slush; or it may continue in well- ‘defined channels. - In a mountainous § 
area of snow cover all of these phenomena may occur simultaneously during 
winter andearly spring, 

- As in other mediums, water will move through the snowpack in a path ot 
least resistance . The rate of movement and the storage of water in snow 
depends on the topography, the temperature, depth, density and the crystaline 
structure of the snowpack. Regardless of its path, before water can pene- 
trate through a layer | of the snowpack, the path in that layer is brought to 0°C 
by the heat of fusion of the freezing water; each snow particle in the path is | 
coated with water, and the capillary spaces filled. Only the excess continues 
downward. Upon reaching an impervious ice plane, the water flows over it 
until a weakness in the ice plane allows part or all to pass through to the next 
‘snow layer below. Cells of cold snow between ice planes may still remain in 
the pack when water reaches the ground level. - However, continued rain and 
melt will eventually disintegrate all ice planes, change the crystaline struc- | 

ture a and render the entire pack isothermal and saturated to capacity. At this 


stage, the water will move downward almost aad snowpack | wil 
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S when Wo is the initial snowpack water equivalent in inches and Tg, is the . 
; average snow temperature in degrees below 0°C. _ If a known amount of water , 


rs. » mental results is considered satisfactory. In general thei: the liquid water 


‘SNOW “HYDROLOGY 


‘yield all input to runoff with little delay. Water freezing within the pack and 
water held against gravity (not including any possible slush water in level 
areas) constitute the retention storage or the initial loss to direct runoff. 
Such water is retrieved only as the snow melts. The magnitude can be com-_ 
puted, and the time delay determined from the input intensity. 

_ The cold content is the heat required to raise the temperature of a “cold” q 

" snowpack to 0°c. Its water equivalent | is the amount of water that would sup = 

‘ply the required heat (80 calories per gram) by freezing within the pack. 7) 

th the specific heat of snow at 0.5, cold content (h) in in is: 


and its water equivalent (W,) in inches is: 


of ( 0° C entered a “cold” snowpack of uniform texture and density, the gen- 
erated runoff appearing at the bottom of the pack would be less than the inflow 
by an amount equal to the water equivalent of the cold content (Wo) plus the 
hygroscopic and capillary capacity of the snow 
_ ‘Like ‘soil, the snowpack is capable of supporting hygroscopic. and capillary 
water in amount varying with the mass of the pack, which is afunctionof ~ 
size, ‘Spacing and shape of the snow crystals. The density, as an integrating 
factor of the structure of the snowpack, appears to be an index of the amount > 

| of such liquid water. Laboratory experiments with melting snowpack of 0.43 
tol 0. 56 g/cc, after complete drainage, showed retention of liquid water ul 


studies are known to show the order of magnitude of the free water holding — 
capacity of the winter snowpack which has generally an average density less 
than 0.43 g/cc. The light winter snowpack with small grains has greater s sur- 
face area to be wetted but larger voids which may not be capable of supporting 
capillary water. ‘Light dry snow, of course, can not remain light after drain-— 
age. When rain falls on it, the snow immediately settles and reaches a den- 
sity near 0.25 g/cc after which the settlement rate is relatively slow. Since” 


snowpack of relatively light density, the lower limit (3 percent) of the experi- : 


in which S¢ is the water stored and f''|, is the liquid-water- -holding capacity 
of the snowpack in 1 percent of W. In the winter, a value of of {Mp = = 0.03 should _ 


prove satisfactory, and in the 1e spring for the daily thawing of the night crust — 
a value of = = 0.04 or 0.05 is 
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in n which D is 7 ain and V is the velocity of water in the liaise The 
velocity of water ina ripe | Snowpack is: relatively large compared | to t the ~~ 
velocity in a cold and dry winter snow. Ina cold snowpack the time of travel 
(t) depends on the intensity of the inflow or rainfall (r) on melt (m), the re- 
tention storage capacity, and depth of the pack. ‘That is, i at 


Where S; is the a amount _ shown in equation (34), and t is the total time elapsed 

Be. from the beginning of rain to the time of appearance of appreciable runoff at 

- the ‘| bottom of the pack. 1 The relative magnitude e of tl the retention storage, (Sy) § 
time delay (t) are shown by the following ‘examples: 


- (1) Given: a cold winter snowpack at 7 ,000 ft. elevation where the snow 


depth is 50 inches with an average density of 0.30 g/cc and temperature of | 

_-5° es rainfall and melt of 2.40 inches is forecast for the 24-hour period at 

an average rate 0.1 inch per hour. ‘Find Wo; the water 

- i ‘Using tp - = 0.03 and V =‘ 40 ins. per hour, & = = 0.00625 (15) ) 5 +0.03 / (15) + 
0. 00625 (18) 57» 0. 93 ins., leaving 1, 47 inches s for runoff which will — 

to appear the at this. after t = 0.10 *40 10. 5 hours from the 


time of of water under the pack, 0: in 


pe “70. 0, 10)) is in the snowpack, moving downward. If this was tne i Rallies rain 
the snowpack had experienced, t>10. 5 because of its detention storage ¢ capaci 


ty in the form of several inches of slush at the bottom of the pack 2 and perhaps 


(2) Given: an open site. ona clear s spring night, the top 6 inches of a 
- snowpack cooled progressively to an average temperature of -3°C by 0500 
hours (sunrise); this cold layer and the remainder of the pack had a density 
of 0.50 g/cc; below this crust the pack was drained of water except for its 7 
retention storage. From solar radiation and temperature forecast, it is esti- 
mated that the melting will begin at 0530 hours, its intensity varying with time, 


between 0530 and 1100 hours, as follows: m =0.01t t2 in/hr. Find | 8p. _ Using 


oa -0. 05, the r night c1 crust will require Sy = 0. 0.00625 (6) (0.5) + 0. 04 (6) (0.5) e7 
0. -18 in. It can be ‘seen | that ripening the snow crust each day requires an ap- 


-preciable < amount. of the daily melt. 


Note that the average temperature and the water equivalent of the snow- 


pack must be known or , before retention storage and time delay 
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HYDROLOGY 


Application of Principles within Columbia 

_ The principles ont procedures developed as a result of the Snow mennineate: 

-is |) tions have assisted in solving problems in snow hydrology in those areas of 
ren, the mountainous western states investigated by the Corps of Engineers. oa 


are not yet completely ‘solved. In the Columbia Basin, use of these principles 
a8, has assisted in obtaining seasonal volume runoff forecast procedures, stand- — 
ard project floods, spillway design floods, and snowmelt relationships and in- 


q _ The Detroit project, located on a tributary to the Willamette River, is a ae 


‘multi-purpose project in which runoff from snowmelt floods is detained and 
Ss used for the generation of power. The North Santiam above Detroit Reservoir 


_— is located on the ' west slope of the Cascade Mountains about 60 miles south - a 

a east of Portland, Oregon, and has a drainage area of 438 square miles. Ele- z 

off at vations range from 1200 to 10,495 feet above mean sea level, the mean eleva- 

2 (Sr) tion being 3718 feet. A large percentage of the area is comprised of valleys 

+e and ridges with steep ‘slopes, anda heavy stand of coniferous trees covers _ " 

now most of the area. Soil cover is relatively thin, but there is considerable duff 

> of and litter under the heavy forest canopy. - The climate is characterized by wet, 7 

d at 4 moderately cold winters and dry ' warm summers. _ Normal annual precipita- a 

liva- § tion over the area is estimated to be 82 inches, ranging from below 70 near lp 

af 5 the Detroit project to well over 100 at the crest, where snow accumulates to _ 


great depths. About 90 percent of the precipitation occurs during the October 
+ The water balance method, discussed earlier in this paper, was wed for Ss 
sin ; » forecasting the volume of runoff from the basin above this project, and the a 
» procedure used is considered representative of a method adaptable to an 
— ) average project basin. Several of the steps taken in this procedure deviate 
re ) from the procedure derived for a laboratory area. No estimate was made 
for losses by interception, since these are small compared to the total rain- 
fall, which in turn is not accurately known, having been estimated by summing — 
the generated runoff and the calculated evapotranspiration loss. - Likewise, no 
estimate was made of gage- -catch deficiency due to the wind. | yo ee 
_ Precipitation, snowfall, streamflow, air temperature, and snowpack water 
equivalent data are available for varying periods. The only y adequate temper- 
ature record available is that at Detroit, necessitating the use of lapse rates 
to obtain estimated temperatures at higher levels. Precipitation data are _ 
available for five stations, one of which was begun in 1909. Water ‘equivalent — 
is measured at four snow courses having records since 1941. Because of _ 
regulation during the construction of Detroit Dam, the streamflow record | 
subsequent to 1951 was not considered usable for this investigation, and the | 


period of record suitable for study extends from the beginning of snow-course — 


iS esti- 
ith time 
.5) = 4 ‘The snowpack remaining on the basin as of August 31 each year is negligi- 
an ap- § ble, and the reservoir is filled after the threat of a rain-type winter flood has a 

oa largely passed, so a reliable forecast procedure for the runoff available from 


February 1 through August 31 was needed. Forecasts covering the periods 
March- August and April-August were also needed. Equation (4) is used for 


all periods. _ The basic equation used for all is as follows: 


, 1959} ASCE Sem, 

— 

gt 

= 

— 

3 

A 


1 January, 1959 
in which Q is the available runoff the the date o of the 
com to August ; 31. = P is precipitation, W1 and W2 are the initial and final ~ 
‘snowpack water equivalents respectively, and Qo : are the runoff volumes. 
under the initial and terminal recession curves, and L is loss, all values 
being estimated as inches of depth over the drainage area. ‘ The final snow- 
pack water equivalent, Wo is equal to zero in this case. The e snowpack — 


water equivalent for the da date of the forecast (W1) is ‘obtained from a a plot of 
the water- -equivalent at each of three key stations, at known elevations, from 

_ which an indication of the water equivalent for | each 10 percent of the a 

_ area is read. The average snowpack water equivalent is obtained by appli 
cation of a factor derived froma study of the relationship of precipitation, % 
loss and runoff during the 11- -year period of study. Qa is the recession aoe 


Pescara corresponding to the discharge on the date of the forecast and yy is” 


climate, and estimates losses by an empirical formula using only air tem-_ 
perature asa variable, with adjustments for precipitation experienced and 
for the latitude of the area considered, which indicates the length of the oe. 


pratep of the water ‘wale during this period. Values of precipitation may oo” 
-assumed on the basis of normal or critical conditions during the spring period 
| 8 someday it may be possible to utilize long-range precipitation forecasts. 
a ) On February 1, the initial snowpack water equivalent Wy amounts to only 
aon 25 percent of the total of of Wy plus P; this forecast is is greatly dependent — 
_ upon future climatological events and is generally to be considered a prelimi- 
- Mary estimate. On April 1, the value of W is about 60 percent of the total of 
_ Wy plus P, and the forecast of this date is reliable enough to be of major as- 
sistance in the operation and refilling of the Detroit reservoir. 
_ The following is the Perens 1- om 31, 1958 runoff forecast of the 


water balance ‘method: 


Recession runoff volume for observed 1 


= runoff volume for normal August 


Feb 1-August 31 off fo t, 
Fe e bruary 1-Augus runc forecast, 


—" ony Dam, now under construction on a tributary to the e Willamette by 
River, is to be a multi-purpose project storing water for flood control and 


j “to this, in order to determine the amount of flood control storage that snout 


be provided if eeateinies feasible. This flood is one which would be ie 
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a ‘exceeded only upon rare occasions caused by a onteenaen falling upon a snow 


covered basin. It requires the combination of rather severe antecedent con- 


= and contribute to the flood. The largest flood recorded at the dam site oc- ror 
ace | curred in 1953, and this rain-on-snow event was considered to precede the 
nes 
is a The 210 square - -mile drainage area above the Cougar dam site, elevation 
a ‘qhes feet ean sea level was divided into five elevation bands of about 1000 
| feet each. Analysis of snow surveys and the climate of the region indicated = ¢ 
t of 4 “that the basin could reasonably be expected to have a 100 percent snow cover | 
from with 8.2 inches of water content at the dam site and an increase of 7.5 inches — 
sin — of water content for each of the five bands. This is considerably more than a ; 
oli- [} the normal and could accumulate only under temperatures less than normal. “a 
; It was conservatively assumed that the basin snow cover was at 32°F and re- — 
_ quired only hygroscopic and capillary water of about 2 percent ‘by weight be- _ 
he is fore yielding runoff, amounting to an initial loss of 0.55 inches, part of which 
Was regained as the snow melted. Melt due to short wave radiation and =a 
ground heat is insignificant during a rainstorm and was omitted. The snow - 7 
melt was computed by zones, using a melt rate of 0.08 inches per degree- -day 
above 32°F applied to appropriate air temperatures for each zone. The melt a 
rate conforms to that previously determined at the Willamette Basin Snow Ba, : 
Laboratory for the condition of rain- “on - “enw in heavily forested areas. _ 


} and a total volume of 194,000 acre- feet. § Storage of 154 ,000 acre- feet of run 
5 off would be required to control this flood, and a usable capacity of 185 ia 
acre- -feet is being 
Some of the variables for clear- -sky snowmelt are neglected in this compu- a 
tation. The heavy forest on this basin markedly reduces the wind speed varia- 
tion over the snow surface. Therefore, a wind variable is unnecessary, and _ 
evaluation of convective melt can be accomplished by use of air temperatures | 7 
alone. _ With intense rainfall the air would be saturated, so that the dewpoint 
temperature would be equal to the air temperature. For this reason, conden- 
sation melt may also be evaluated from an air temperature | function alone. aed 
Rainfall over the basin during the most intense 6-hour period was estimated 
to be 3.26 inches, and snow melt during this same period was 0.72 inches. _ 
Total rainfall for the design storm amounted to 10.75 inches, and snowmelt _ 
totalled 3.64 inches. The temperature at the 3900 foot level was assumed to > 


2. Several spillway design floods have been computed for projects in the Col- 
umbia River Basin using the procedures and relationships that have been de- _ 
veloped as a result of the Snow Investigations. Such a flood was developed for , 
the Libby project, on the Kootenai River in Palade and just recently for | “— 
projects on the Salmon River in Idaho. These determinations require a con 
siderable amount of detailed work, and they will not be described here. Bas 
- A very - important application of snow hydrology in 1. the Columbia Basin is 
evolving in the day-by-day use of data gathered by electronic means which 
must be interpreted and used in producing a forecast of flood flows and gage 
heights in the ensuing short-term period. Reservoir regulation for flood a 


nilar control and flood fighting operations at local critical areas and in the reach 3 


should of leveed | river downstream f, from Bonneville Dam \ will be assisted thereby. 
Because a forecast must be completed within a short time after the raw 
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“become available to be of use, use, the hand or "| ‘machine computations enter- 
ing into them have been restricted to a relatively few. With the advent of a a 
_ high speed digital computer, and a routine written for the Columbia Basin, a 
7 - forecasts will be made for a large number of small tributary basins, based © 
upon antecedent weather and forecasts of temperature and precipitation to be 
expected in the next several days. . The same routine will route the flows < 
_ = the small tributaries through lakes and channel storage and will furnish 
plotted hydrographs at key points. ‘Such a routine requires the derivation of 
"indexes and coefficients pertaining to each of the small tributary basins. se 
a a One of the tributary basins selected for analysis was that of the Boise | 
_ River above the Twin Springs gaging station, where the 1954 and 1955 spring” 
snowmelt hydrographs were used. Boise River above Twin Springs, 
Idaho drains an area of 830 square miles of the Sawtooth Mountains. “The | 
area has deep valleys, steep slopes and narrow sharp-topped ridges, ranging 
in elevation from about 3000 to 9000 feet above mean sea level. "The stream — 
system and instrumentation are shown in Fig. 3. Forest cover consists © i, 
- principally of conifers which are estimated to have a Shading effect on 30 per- 
cent of the basin area. _ Snowmelts for each day of the periods were ‘computed 
by the thermal budget index and the flows resulting from this released water > 
‘plus the precipitation for that day were routed to the gaging station for com-— 
parison with the observed hydrograph. This comparison for the May- -June 
_ 1955 period is also shown in Fig. 3. _ The index derived for the snowmelt over 
_ where M is the daily snowmelt in inches over the basin, Tmax is the daily 4 
- maximum 1 temperature, at Idaho City, in degrees F, and G is a radiation _ 
; parameter, in langleys (equal to absorbed shortwave radiation plus es esti- — 


m It will be noted that one of the terms in the above equation requires that — 


a forecast of the radiation to be experienced on ensuing days must be availa- } 


ble. Radiation is measured daily at Boise, Idaho and these data were used in fj 


deriving the equation, but large parts of the Columbia Basin do not have such | 
‘measurements available. | _ This refinement may be used sometime in the fu- 
ture, when more reliable forecasts of the temperature and precipitation and» 
other weather factors are possible and when data on radiation become more 
widespread. _ For the present, a simple maximum daily temperature index of 
‘snowmelt is being used, which provides fairly reliable estimates of snowmelt 


runoff for day-to-day forecasts of streamflow. 


__ Another application of indexes to the predication of flood volumes resulted 

j . The method is based primarily | upon the rela-_ 
tionship between winter runoff of low-elevation drainage basins inwestern — 

_ Washington and Oregon, and the spring snowmelt runoff of the Columbia | 
River at The Dalles. Indexes of winter temperature and spring precipitation — 
4 are included in the forecast procedures as secondary parameters. The use — 
of low-elevation winter flow as an index is confined to regions where the low- 
elevation and the high - -elevation areas have a common source of moisture. ae 5 

4 Such a situation exists in the region composed of the Columbia Basin and © e 
western Washington and Oregon, the entire region being well centered in the 


belt of prevailing westerlies. Moisture is carried i in a generally | eastward 


direction from the Pacific Ocean, and the amount deposited is largely 
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air, reflected 1 by both winter ‘streamflow from low elevation areas and the . 
a accumulation ¢ of snow at high areas. Several tests indicated that a correla- 
tion of winter flow from two streams each in Washington and Oregon would 


yield results as accurate with regard t to historical data as those which use q 


- dependent upon the rate of flow and the pasclatabiie moisture content of the a 


» precipitation and ‘snow course data for the pr incipal indexes. This 1 method — 7 
a does not produce any information as to the distribution of the snow pack over 

the smaller interior | basins « of the Columbia River above ' The Dalles, Oregon, 
but it does produce : an acceptable forecast for the basin as a whole quite cane. 


n th lati d. 


SUMMARY AND CONCLUSIONS 


‘The : significant developments in applied snow hydrology which have re- 


sulted from the work of Snow Investigations research and analysis are: a 


A point snowmelt equation based on heat supply from: mab and longwave 
"radiation, convection and condensation. 
_A general snowmelt formula applicable to any eileen area when over e 


coefficients characteristic of that basin are evaluated for the degree of ex 
q posure to wind and radiation. vine 


Relative magnitudes of melt different sky and forest 

Numerical appraisal of the storage and delay effect of the snow pack on 


Runoff forecast on -balance and on a coastal 


In addition to the numerous investigations were 


ae which were of a more basic nature. _ These deal with the theoreti- : 


cal aspects of snow physics, precipitation, snow cover, and soil moisture _ 
measurements, forest effects on snow accumulation and evapotranspiration 
loss, heat exchange | between the ‘snowpack and its environment, the distribu- 
tion of the snowpack in mountainous areas, and the water balance in areas ( of 
Our knowledge of snowmelt is not complete, but we know enough now ee 
supplement the degree-day approach by including humidity, wind speed, and > 
radiation values even if they must be partly estimated, particularly | for _ 
_weather melt in open and partly forested areas. Only ‘during a rain in .. 
dense forest, where the wind speed and the solar radiation are small and 
nearly constant, can the temperature alone define the snowmelt properly 
Direct application of these principles is now being made in estimating de- 
_ Sign floods for determining the capacity of flood control storage projects and 
river channels, and the related problems of conduit capacity, revetment de - 
sign, levee heights, , and day - -to- ~day operation of the reservoirs . Seasonal 


ae ~ The network of stations is being expanded to include more eae in 

formation and better areal coverage. Radiation data are observed at 
few stations and this class - observations should be extended to cover more 
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area. snow fields. in the | high are and observa- 
tions of climatoligical events in this large and important area are difficult to 
obtain but steady advances are e being | made by use of radio-reporting facili- 
ties. The distribution of snow cover during the melt period is also a ob- 
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THROUGH CIRCULAR WEIRS* 


Closure by. é..C. 


_J. C. STEVENS, 1! 1 M. ASCE.—The is gratified th six engineers 
well versed in ‘mathematics 'S and hydraulics: took time to send i ‘in discussions 

’ The experiments by Mr. Manabar are very interesting even though his 

ra pipe > was ‘small. _ The ‘sharp edge of f the inlet: of a pipe on a steep slope is not 

® the same as a sharp- edged circular weir through which free flow obtains. 

© When partly full the inflow was partially submerged and this is correctly as- 

© signed by him )as the reason his average flow coefficient of .52 was less than 
the author’s average of .59. . His Fig. 2 indicates that the stage- -discharge re-_ 
lation is virtually a straight line for higher values of H/D. This is in con- a 
formity with the graphs for the laboratory flows listed in Appendix A. a. 

» The head on circular weirs need not be limited at all to HCD for the weir 


may be surcharged; also it may be submerged and it is often necessary to 


calculate the flow under these conditions. 


_ For surcharged ‘circular weirs there are two conditions to be met. In the 
range H<D flow is readily calculated from the elliptic formula with the aid 

of the tables in the author’s paper. _ When HDD the weir becomes an orifice — 
for which the common or orifice formula is s applicable, the head being measured 
from the center. Fig. 3 shows the calibration of a 6- inch sharp- edged circu- 
lar weir by E. R. Dodge to a head of H=D (Appendix A). _ However, , his thesis(3) 
gives the flow for | ‘surcharged heads. ‘The weir was in the outlet end of and 
concentric with a 15-inch pipe in which the approach velocity was negligible. 7 
The av erage flow coefficient of -607 given in Appendix A for this weir was : 

u used with the Equation (3) to define e the solid line Fig. 3 (a) for H>D. For 


surcharge heads the same coefficient was used in the orifice formula viz 


= AV2 ‘gh A =.197 sq ft and h is the head on on the center of the 


given the effect of surcharging the 1-foot weir ‘from 
I omori’ s ‘experiments in his Table B. Fig. 4 is a graph of the flow through — 
I the Omori weir calculated by Eq (3) up to H/D = 1.00 using the average flow 

} coefficients of the last column in Table A. For the surcharged head the cor- 
responding a average coefficient for the w weir, .596, , was also used. Note ethe 

+ smoothness of the transition of flow weir to orifice in both examples which | 
are well substantiated by laboratory experiments. 
For submerged flow the Villemonte formula (4) may be used. For a we weir _ 
of any kind subject to submersion in a stable channel observations will soon At 


a. Paper 1455, December, 1957, by J. Stevens. 


wd Cons. Engr., Stevens & Thompson Engrs., » and Leupold wai Stevens 
ot Instruments, Inc., Portland, Ore. 
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ms “nel with a proportional sharp edged weir at its outlet. Baan 
The general form of the Villemonte Formula is 


‘Where Q is the observed flor flow -cis. 
the flow corresponding to the headwater (above) 


Qp, the flow corresponding to the tailwater 


aconstant 


‘an ‘exponent to be determined 
‘Both Qq and Qb are taken fi from the free flow curve of Fig .3 (a). a For most my n 


7 -sharp- -edged weirs k is practically 1.0 and m has the average v: value of .385 
_ found by Dr. Villemonte’s laboratory experiments. A nomogram is readily 
i _ constructed for the Villemonte formula for any set of stable conditions. (15) 
i In Fig. 3 the dashed line shows ; the effect of submerging the weir- -orifice. 
: This well illustrates the versatility of the circular weir-orifice. What other 
id of open channel measuring device lends itself to reasonable accuracy 
_ throughout a head range from zero up to a hundred or more feet? Even the 
Venturi tube is not accurate for low rates of flow. 
_ The author is grateful indeed for the errata in Mr. Blaisdell’ '$ discussion. 
_ As one grows older his facility for errors increases markedly, especially " 
4 with the slide rule. . The flow experiments have been checked and revised — 
i: tables of Appendix A, where errors were » found, are appended hereto. * ie 
a table of Appendix C is revised in its entirety in order to eliminate the 
minor errors in the 4th decimal place for which Mr. Moody takes the author 
a to task and quite justifiably so. . Any table should be correct to the last deci- 
mal. Mr. Moody’s admirable table of values of Ato 4 decimals with consis- 


— 
oOo 
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tent 2nd is a step in the right direction but interpolation may be necessary. > 
The author believes that ‘ke ‘should be carried to 3 decimal places. tia: 
Mr. Kolupaila’ s table gives the theoretical flow directly for a 1- foot diame- he p) 
ter weir. The argument x has been substituted for k2 in the author’s Eq (3). - 
The theoretical flow through a circular weir of any other diameter D is ob- , 
tained by multiplying the tabular value by 0 


answer to Mr. Blaisdell and doubtless others who noticed the lack « of in- Pa 
formation regarding the approach conditions for the weir data in Appendix A, whic 
references for the experiments now include this information. = Ff | 

The author knew nothing about the Byrd and Freidman Handbook of Elliptic r 
Rimi for Engineers and Physicists, Moody’ Ref (1). He has since pur- 
chased a copy and was quite non-plussed ” find included 2 errata pages of 7 
fine print — - about a hundred corrections. In the sweet -bye-and- ~bye perhaps 


_there’ll be no errors - “Tis a consumation devoutly to be wished. er ae 
an In weir computation the orthodox method followed throughout the entir is | 


them with observed discharges and correct computed flows by a factor 
which might be called a coefficient of ignorance. esses 
For all sharp edged weirs that coefficient is some > measure of the effect of 
the sharp | edges distorting the stream lines as. water passes the weir. Mr. 
Blaisdell’s use of his formula (a) for circular weirs is unorthodox but cer- 


tainly utterly one even measure H from 


Gefine a relation between headwater and tailwater. To illustr 
J 4 Oe shows such a relationship (similar to one found for the outlet eg 
ven | 
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any ¢ arbitrary detute, reduce or increase a by a a constant or a percentage, 
even use the value of g from Mars or Venus and obtain “goofy” discharges 7 
(Qnuts) computed from (A). Parallel observed discharges would then deter- 

| mine a series of coefficients that would define a “goofy” curve but from which 
the true discharges Q could doubtless be 
= The average coefficient C= .365 (Fig. A), he admits gets top-heavy beyond - 
H/D = = .8. Dividing the scattered values of his c by the author’s correspond- — 
ing scattered values of C (Fig. 2) he obtains the tantalizing curve of Fig. B _ 
which the author at first mistook for an improved set of discharge coefficients _ 
applicable to circular weirs. However, Mr. Blaisdell, by letter, disclaims 


any such intention and states “the values in my Fig. Barea measure of the 


incorrectness of my approach, part of which I suggested might result from 
my neglect of what Professor Mavis calls ‘shape factor.’ I worked it out and 
compared my results with yours as a matter of ‘interest. . . I thought others 

The data on discharge coefficients by. Messrs. . McPherson and 
‘Kolupaila are shown in Fig. 5. Curve I indicates the average coefficients of 
the 4 e experimenters given in Table A of Mr. McPherson’s discussion viz 


Ormori, Staus, Staus & Sander, andGreve. + 
Curve II gives Dr. Kolupaila’s ——— formula’ which he > slightly modi- 


fied from | that of Staus | as follows: en 
0. 570 + 0. 


100x 
where x = _H/D. The Staus formula used 0. 555 in place of 0.570. 
The author’s constant coefficient of .59 fits Curve I fairly well for moder - 

| ate to high values of H/D. The deviations are “9% at H/D = .10, -5% at | -20, 
and practically zero above 


Moody has offered a new and valuable approach to orifice flow by his: 
formulas (3) and (4). In order to avoid confusion it would be preferable to use 
(D/H)2 rather th ‘For theoretical flow through vertical circular orifices 


tie 


Where . Ais the elliptic in Appendix C using — for the argu- 
ment instead of k2 = = H/D, while H is measured from the invert of the fone 
as in the weir instead of from its center. Obviously D/H< 1.00. 
“Fig. 3 shows by a dotted line the results of using (21) for the flow through 
Dodge’s 6-inch orifice, using the average weir coefficient of 607. (Appendix — ; 
A). Fig. 4 shows the corresponding flow for Omori’s 1-foot orifice in Table BO ; 
using the average coefficient for Omori’s weir experiments of .583 given in 7 
the first column of Table A. It appears therefore that the average poe 
coefficient found for the circular weir may be used as the average coefficient — 
for the same weir when it is surcharged becoming an orifice as is — , A 


by the four examples herein given. More dunes experiments on sur- 
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iL. Cone, V. M. Circular Weirs. Journal of Agricultural Research USDA 4 


March 6, 1916. The approach flume was 20 feet long 10 feet wide and 6 


center of the weir. 
, F. V., Flow of Water through Circular Parabolic ann: st Sepia 


"Vertical Notch Weirs, Engineering Bulletin, Research Series No. 40. pee 
Purdue University, Lafayette, Indiana. _ This bulletin also contains a 
bibliography with summaries but no mention of European a 
_ The approach channel was 16 feet long and 5 feet ' wide. The depth below 4 
the weir is not given but it scales 2-1/2 feet below the invert of the largest 


3. Dodge, Eldon R. A Thesis Submitted for a Masters Degree, University | of 
Wisconsin, Madison, Wisconsin. 1935. (Not published). The sharp- 
edged circular weirs reported were in the end of two 
sizes of tile sewer pipe as follows: : 
Weir Nominal D Pipe I D nal Pipe I D : 
1/2 8-inch inch 15-inch: 


‘5. Thijsse, J. Th. 2. Meters Director Hydraulic Laboratory 
_ (Waterloopkundig Laboratorium) Delft, The Netherlands. Three sizes of 
- flumes were used. Weirs were not rated beyond H/D = .75. For the ‘ 
_ 2.36-foot (72 cm) weir the nominal dimensions of the flume were, width 3 
feet, depth below invert of weir 1.0 ft length of 9 ft. There is strong evi-— 
_ dence that the increase in the coefficient for all 3 sizes is due to velocity 
approach. However weirs were calibrated in place. 


- For the 1.77 foot (54 em) + weir the flume dimensions were width 1 2.2 2 feet, 
depth below invert of the weir feet, length 7.0 feet. 
= “For the 1.181 foot (36 cm) | weir the flume length was 4.6 feet, width 1.5 


14. Submerged-weir Discharge Studies by James 


News- Record Dec. . 25, 1947 ‘and 


af 


6 add the following references: 


‘Proportional Weirs for Sedimentation Tanks, by 3.4 Proc. 
"Appendix A - Substitute the following tables for those of the corresponding 
Author an¢ and diameter. 
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Same Pp 24 Eq wes 
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Page 5 line substitute: 


The by V 


. M. Cone are near the an line 


made his experiments in 1916 at Fort Collins. mm & 3 
7 gives A to 4 decimal places for values of _ — a 
a 
es 


January, 1959 


— 
0.0000 0.0000 0.0000 0.0000 0.0000 2.000 9.000 0.0001 
0.0003 0.0004 0.0004 0.0005 0.0006 - 0.0007 0.0008 0.0008 
0.0013 0.0014 0.0015 0.0017 0.0018 0.0020 0.0021 0.0023 
+0026 0030 0.0032 0.0034 0. 0036 0.0038 0.0040 0. 0042 0.0044 


0073 0076 0.0079 0. 0. 0. 0088 0.0091 0. 0094 0.0098 0. lo. 
0.0104 0.0108 0.0111 0.0115 0.0119 0.0122 0.0126 0.0130 0.0134 
0.0146 0.0150 0.0154 0.0158 0.0163 0.0167 0.0171 0.0176 0.0180 
0.0185 0.0189 0. 0194 0.0199 0. 0203 0.0208 0.0213 0. 0223 
-09 0.0233 0.0238 0.0243 0.0249 0.0254 0. 0259 0. 0265 0.0270 0.0276 0. 0261 
+0287 0.0293 0.0298 0.0304 0.0310 ons 0. 0322 0.0328 0. 0. 
0. 0.0353 0. +0359 0.0365 0. 0378 0.0385 0.0391 0.0398 0.0404 
0.0411 0.0418 O425 0.0432 0.0439 0.0446 0.0453 0.0460 0.0467 0.087% 
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‘Discussion by C. O. Clark 


cs CLARK, ee M. | ASCE.—The mechanical analogs of aid in graphical 


flood ‘routing ‘presented are are of interest to all who must manage the natural _ 


waters of rivers and streams, either by scheduling and alloting the time of 
_ availability o} of beneficial use of that water which is harnessed and under con- a 


trol, or in warning of the unavoidable rampages | of the public enemy when it 


‘The philosophy by which the author searches for his mechanical analogs, - 
_ finds the dampening effect of storage to be kin to the effect on the move- 
. ment of a rear auto wheel of the distance between the wheels, makes a 
ing reading. ~ Such refreshment is needed by those who find the need for think- 
ing and the opportunity for it completely eliminated by the codes, rules and oe 
: regulations, and “accepte ” formulas which make up much of the necessary _ 
pieces of routine engineering work. There are excellent examples of apple 
logic and scientific observation : in the manner by which each of these analogs 
has been developed and built into a solid tool. There is food for thought in the 7 
fact that the author who conceived and built electrical analog computers for 
the solution of these problems, has now turned to mechanical analogs. ie 
we One of the analogs is as small and as subtle as a sliderule, and is as per- 
sonally possessable a asa straight edge/or a pencil. it is a truly Euclidian y 
— tool . Even more, by an as yet unproved modification of this tool, the. author 
_has opened a promising door toward a variable, non-constant, kind of solution | 
3 ofa problem which exists in the nature of streams as they : are, but it which — a 
7 "solution has defied the limitations of brains, two- dimensional | paper and the cae 
7 timeclock up to now, namely the class of problems involving variable relations oll 
between storage and discharge, variable > lag and storage. 
la The writer is interested in the template type analogs, but this does not de- 
tract from the more mechanical ones, which are certainly ingenious, probeb- 
| ly accurate, and based on sound analogy. . This interest is very strongly in- 
, fluenced by the fact that the underlying principle of the templates, which oa 
_ Kohler now applies to long reaches of the main stem of a river system, is the 
same principle which underlies an application to an entire watershed to de- 
velop the type of flood hydrographs which the writer could call unit hydro- 
graphs, if they had a unit duration of rainfall, but must call instantaneous 
hydrographs | when that unit is zero. This essential principle is that two steps 
_ of simple analysis: are first found(1) to be essentially comparable to the more 4 
complex and simultaneous action of aquatic storage as it occurs in nature vand 
as McCarthy approximated it in flood routing formulas developed for the _ 


a. Proc. Paper 1585, April, 1958, by Max A. Kohler. 
Hydr. Engr., Tulsa, Okla. 
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author calls it lag routing, or r lagging and then routing storage mome 


; which he explains is bodily translation in time of a hydrograph without stg and d 
ce tion, followed by a a distortion or attenuation much in the manner effected by routit 
storage in a large level reservoir where storage is strictly a function of in- 4 = 
_ stantaneous outflow. | This is a mental analog, of course, because the one kind a ] 
apply 
ofa waterway in which none of these mechanical analogs is needed is the level - effect 
reservoir of great depth. Nevertheless, it is a correct analogy andit is the 
Limiting case of a class of problems which vary from the headwater problems ie 
 gtrea 
of sewer inlets where translatory effects have been c considered large and pond- tion 
ing effects negligible, down through the streams which embody both in every — a expel 
_ Proportion, onward to the large > reservoirs, great lakes, and then to the ace 


‘seas in which the largest w: waves made by floods are but Tipples on a surface 
_ translated in negligible time and at velocities compared 1 with that of sound, | 
and involving: no appreciable lateral movement of mass. This idea i is not pren 
‘Translation or lagging was built into flood calculation and thought in the sO- 


called method of decades ago, but not originated by Ramser_ 
comk 
_ problems of drainage at the time that Yarnell and Ramser were young men in . 7 2) 
these fields. | ‘It is as natural fora headwater engineer t to think in "mate 
quen 
Ir 
whic 
voirs, but not so effective in attenuation. . The writer has tried to knit lates 
these upstream ideas of time required for concentration and the downstream © ; upst 
_ideas of storage which must be filled before flow can become steady, into the Z This 
concept of an instantaneous hydrograph(1) which is the least size of hydro- _ McC 
graph that could have a peak and the minimum duration if all elements of the wie 
_ watershed were equally involved. | Hence it is dealing with a minimum volume | upa 
associated with a flow anda _ duration and an influence of each on the other and pone 
7 of upstream and downstream elements on each and both, or the minimum Ze buil 
_ wave which is stable enough in form to establish essentially stead steady flo flow an ant tive 
transmit. it along a sloping channel. into 
_ What is definitely new in this work of Kohler’ s is the correlation of transla- ‘Sent 
tory storage with upstream storage, or lag, the downstream influences with | - suit 
"change in shape | as distinct from that of pure translation. It has been done ; ad 
before, and justified empirically or non-dimensionally in several trial-and-_ \ 
- error or convenient solutions. - But in words Kohler seems to be the first who. Sgr 
has said lag is a function of inflow and attenuation is a function of outflow and | den 
who, in dimensions, has put one on one side of a line measured one direction — = 
and the other on the other side | of the lin line measured the other direction a as bet 
 hehasinhis template offigure7. 4 | 
__ The writer was not prepared to accept this at face value as rational. 2&B But § mil 
man already committed to an empirical usage of that principle, and having sal tive 
the means of experimentation at his disposal, he caused two flood waves to be ; rat 
_ generated by purely upstream forces, within the banks of a known stream, and _ tion 
the lateral inflow which confuses natural flood records. was 
| 
case with most rain-generated waves. The stream is steep enough so that : | a 
aage- -discharge relations are as constant as they can be in ordinary — ~* | 
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— 
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streamflow hydraulics, and the friction factors are so much larger than 
“momentum forces and velocity heads that influence of these things on stages” 
* and discharge is probably much less than it is in the streams for which flood 
routing techniques would be applied. The results are summarized here as) 
‘support for the author’s view (that the Muskingum routing sequations 
apply, i.e. constant time elements and linear storage-discharge relations, the 
_ effects of of flood routing can be suitably approximated by transposition asa 
function of inflow and storage attenuation as a function of outflow) (The 
: stream does not have the form or dimensions which would lead to intensifica-_ 
- tion of wave form, , $0 these forms are neither | confirmed nor denied by this — 
oO premises are common to the selection of the experimental wave forms, 
the author’s lag routing templates, and the writer’s concept of hydrograph “a 
7 analysis and synthesis by basin shaped instantaneous unit- -graphs. a These two a, 
- premises are that, in a stream where time elements of flow are approximately 


a 1) Storage rage effects may be approximated by two kinds of storage in | suitable 
ie which tend to attenuate flow and which do not, and, ye co 


of flow and increasing and decreasing flows | may be ‘approxi - 
wi by small instantaneous impulses of positive or negative change, and 
hydrograph no matter how complex c can be sgpronimates some se- 


a In the writer’s view, premise (1) is ouae of first as amounts of storage — 
- which translates and attenuates, or reductive storage, ‘and some which trans- 
rf lates and accentuates peaks, with the latter apparently a function of inflow or 
2s factors” and the former a function of outflow or r downstream factors. Z 


wise be combined i in combinations of two prisms so that tendancies to build . t 


build up which is possible in some cases and would d apparently valve nega- 
tive values of the K-item). ‘This prism combination is Mr. Kohler’s division — 
into lag-and-storage values. (The writer would caution that both items repre- | 
sent storage, however, and that translatory storage and reductive storage = ; 
‘suitable names which retain the concept of physical storage space required in 
channel to produce the effects involved). 
With the above premises in mind, the experiments here reported were de- 4 
a ‘ean to 1 release two rectangular wave forms of different magnitude by sud- 
denly opening and suddenly closing gates with suitably sustained periods be- 
tween each change to permit flow to become completely stable and steady nm 
between changes and steady throughout the zone of the experiment. ee 
“le resulting wave was observed as river stage at two points, 0.3 and 6.7 
_ miles downstream and at five-minute intervals. Using the two stages of posi- 
_ tive steady flow observed and the initial stage as zero flow, an approximate _ 
rating curve was prepared for each observation point and the stage observa -— 
7 tion converted to approximate steady flow in much the manner which is conven- 
q tional in stream gaging on the steep streams to which these premises can 
of. The ag an 
The interpretation is shown on figure 1 lag 
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case, or when inflow ceased 


are the conventional ones. In this case, unlike natural hydrographs in which 2 
inflow is never steady except at zero, a determination is also possible on the _ 
_ ascending leg, and the values are shown. 
Aa Within the range of many practical uses, the values rising and falling on a 
‘the associated lags are constant. Whether the small differences in the time 
elements ai as printed on the figure are factual or the faults of methods of ob- A 
servation or deduction is not known, and the opportunity to repent the experi-_ 
_ ments suitably has not since presented itself. 
_ The measure of the whole storage utilized in wave passage is the sum of © 7 
‘the: associated lag and recessional values, successively 258, 298, 310 and 346 q 
minutes. As volumns these are respectively 0.36, 0.41, 0. 43, and 0.48 acre- 
= per cis of flow, or nearly 800 acre-feet in the prism of changing water a 


level in the range of these observations. The writer wishes to emphasize | 
that these are physical volumes, and changes in channel alignment or confine-_ 
_ ment could alter them in proportion to those changes. =~ 
‘This is a stream where the proportion of storage along the channel to flow 
‘in n the channel is about the same throughout the range of in-channel flow rates, 4 
as has been confirmed in backwater computations using channel cross- Fs ts 
‘sections. 2 Hence, this is a stream to which constant time elements apply es-— 
‘sentially within the channel and the range of these observations. The writer 
cautions against too rapid extension of this experience into streams where . 
constant time elements do not apply, and where linear | storage-discharge re- 
lations are not a fact. All of the “black marks” on the record of the McCarthy — 
s flood routing method and on unit hydrograph theory ha have been put there by oi 
people who did not first find these constant time elements or linear relation-— 
: ships to be facts as applicable in their subject streams. Where these rela- 
tionships are factual, the theories work very well. _ However, where t the ap-— 
plicable facts are that the storage - -discharge curves are wide and 
loops, and there is no fixed relationship between stage and flow rate in the 
first place, more basic and more fundamental methods must be used. ‘The 7 
-variable lag -variable reductive storage templates illustrated by Mr. Kohler 3 
offer more promise than anything else yet available, but, like many other 
loop relationships, may not be constant from flood to flood and wholly precise “ 


| 


aa 


except in the flood from which the loop or template is derived. _ Nevertheless, — 
‘Since the streams to which such loops apply have other severe sc | 


_ On the positive side from this cautionary expression, the varia variable templates ; 
are e capable of indicating the earlier arrival of very large flood peaks | and the 
_ almost complete flattening of small floods just beyond the stages of bank « over - 7 : 


flow, which phenomena are the characteristics of ny oa larger streams and _ 
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_ probably of all aggrading, silt- -carrying streams (like the Verdigris River 
_ which Mr. Kohler used as an example for his template of figure 7 and illus- | = 
tration of the presented | idea) ). Because t the > template is adapted to the simplest — 

of graphical methods of solution, and involves no reading, tabulation, replotting 
or reprinting, but produces a hydrograph from a hydrograph, it is the most 
promising presently available answer to some difficult flood-time problems. > 
One of these is the solution for many streams almost simultaneously, as in A 


the emergency of flood - time, and where constantly changing data are n— 


»mergency ‘conditions and can be safely done with the templates c care 
fully prepared under the more leisurely conditions of other times. | a 
aa The template is the newest and most unique of recent flood tools. It i is 
destined to rank for several decades with Sherman’s -unithydrograph as 
basic simplification of essential facts for practical daily 
eo Sense rage and the Unit Hydrograph, by C. O. Clark, Trans. ASCE, oe 110, 

1490, fig. Set seq 
; 2 Engineering Procedure as Applied to Flood Control by Reservoirs, aaa 

~~ Reference to the — Flood Control Project, Engineer | School, 

3. Areas, by C. E. Ramser, Journal of 
‘Agricultural Research, Vol. 34, No. 9, May — ag 


4 Engineering Flood Control, X-156, Engineer School, 
Belvoir, 1946, p. 134 from sources of about 1929. 
L. G. G. Puls). 
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by Leo R. Beard 


LEOR. BEARD.- _This paper Mr. Childs. 
hydrologic aspects of a record- ~breaking flood event and their effects on the 
design of flood-c« control . The writer is ‘particularly interested in 


feels that “flood frequencies for stations es 100 years of record will give * 
} fairly: reliable frequencies up to 50 years. . Beyond this point the reliability 
_ drops rapidly and beyond 100 y years the frequency is fictitious.” However, — 
- estimates of extreme flood frequencies or probabilities are necessary to an a 
adequate economic evaluation of project justification, and proper planning 
therefore requires s that satisfactory methods of extrapolating frequency mn 
curves be developed if at all feasible. The writer has briefly investigated 
this problem along lines explored by Messrs. Russell Morgan, Marshall 
-Iakisch and others of the Philadelphia District, Corps of Engineers, in Seite 
studies of the Delaware River Basin. It is hoped that the following brief 
; Beet of this investigation will shed some light on the evaluation of fu- - . 
“ frequencies | of extreme hurricane floods such as those of August 1955. an : 
L 2 Experience in computing flood frequencies for the New England area has -_ 
“been particularly discouraging. 2 While in most of the United States a 
' maximum flood events are apparently distributed in reasonable accord with 
the normal probability curve, there has been strong evidence that the New ' 
_ England area experiences an unusually large number of major floods (that is a 
that the frequency | curves are skewed in the positive direction). Until the — : 
7 1955 floods occurred, there was some debate as to whether the observed ex-— 
cess of large floods has been accidental (due to chance) or r whether it derives 


rest of the 
; ae The author suggests that “ “large floods and the smaller events are not re-- 
lated and should not be compiled in the same array and analyzed by statistical 
_ methods,” and “analysis of floods should be divided into two categories: « : one — 
_ for moderate annual floods, and the second for the larger floods.” The writer 
believes, , however , that those factors 5 causing a major event can also combine 
= as to cause a small event, and consequently that some small events in 
must be related to the large floods. It is doubtful that there is an abrupt dis- 
continuity in channel conveyance factors such that they are distinctly | differ- 
a However, it is agreed that events which are not related (that is, not caused 
_ by the same e general factors) should not be combined in a ‘Single frequency ar ae 


Tay. ‘Since most of the New England floods are caused bye extratropical aes 
a Proc. Paper 1663, June, 1958, by Elliot Childe. 
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"January, 1959 
= saci but a fom are caused by tropical cyclones (hurricanes), it has been a 
suggested that events be segregated into two groups on this basis. Inthis  _ 
: -= perhaps theoretical distribution functions used successfully in other re- 
gions might be _ applied in the New England area. .. Such functions are > requisite — 
= _ to accurate evaluation of extreme flood frequencies. 
be _ _ This idea which was developed in connection with Delaware River studies, | 
- was tested by | the writer in a study conducted under one of the Civil Works» a ; 
Investigations Projects of the Corps of Engineers. Twelve long- record 
stream gaging stations distributed throughout New England were selected d for { 
hi 2 study. The period of record considered was limited to the last 40 years, 


because data on hurricanes during earlier years are not reliable. Maximum 
daily flows were selected in three series as follows: 


a Annual maximum non-hurricane floods: 
Annual maximum hurricane floods 
. Annual maximum floods, regardless” of ‘type 


The list of hurricane floods used was compiled from various hurricane re- 


ports. | The various reports available \ were ‘hot consistent, and iti is felt that, 


been listed. Hurricane | floods used in the study are as follows: 


Flood Dates by 10-ye: -year Periods 
25 1926-35 1936-45 —-1946- 55 oa 
23 Jul 1916 Oct 1927 19 Sep 1936 30 Aug 1949 
1 Oct 1920 Aug 1928 21 Sep 1938 Aug 1950 
Oct 1923 Oct 1929 Aug 1939 (12 Sep 1959. 


27 Aug 1924 (17Sep 1932 Oct 1939 Sep 1952 
1 Oct 1024 24 Aug 1933 2 Sep 1940, 15 Aug 
18 Sep 1933 Sep 1940 Sep 1953 
Jun 1934 17 Oct 1943 31 Aug 1954 
9Sep 1934 2 Aug 1944 12 Sep 1954 
5 Sep 15 Sep 1944 16 Oct 1954 


Oct 1944 Aug 


_ Since the frequency curve of hurricane floods is based on a ‘events 
a — the years of record, plotting of the curve was based on the percentage of 
- floods rather than on the percentage of years. In order to convert the fre- — 
a quency in percentage of floods to the frequency per hundred years, the indi a 


cated exceedence frequency | must be multiplied by the ratio of the number of | 


cane frequency curves is illustrated in Figure 1. 


was found that, in the case of hurricane floods, the plotted data fit the com- 
puted log-normal curves reasonably well in most cases, but occasionally sug- 


= an upward curvature of the ‘frequency curve. The plotted data of the non- 


= floods used to the number of years of record. This adjustment of the hurri- | 
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January, 1959 
- order to summarize the data for all onan: stations in such a way as to 


cess of the indicated 100 ~year \ value were > counted. Similarly, the floods that 


forth until the entire range of eqeeney was covered. Comparison of ‘these 

aggregate numbers with those indicated by the straight line frequency cu curve 

is shown in Figure 2. It will be noted in Figure 2 that the hurricane-type - 

floods evidence a slight upward curvature throughout the > range of the fre- | > 

quency curve. The non-hurricane floods evidence a sharp ‘upward curvature 

only in the upper range of floods. Since it was known that the largest non- 4 
_ hurricane general flood that has occurred in New England in the past 300 14 


years - is that of 1936, it was felt that perhaps this one flood may have been 7 
unduly responsible for the upward curvature shown by the non-hurricane _ 
floods. In order to examine the effect of this one flood, the 19. 1935 flood events 
were subtracted from the frequency count and the aggregate points” replotted. 
This comparison, also shown in Figure 2, is is somewhat better, but there still 
a is a tendency toward upward curvature. 
When floods are segregated according to > cause or typ type, ‘and where the 
meee frequency curves of annual maximum events are plotted as in this | 
study, it is often n desirable to obtain a frequency curve of annual floods re- ir 
gardless of type. In order to do this, the following formula — in the 
Philadelphia District, Corps of Engineers, is employed: 
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= per hundred years, annual maximum hurricane 
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frequency per hundred years, s, annual maximum floods, 


Combination of frequencies using this formula is illustrated in Figure 1. i 


_ As might be expected from the experience with the separate hurricane and 
non-hurricane flood frequency ¢ curves described above, combination frequency 
curves “computed in the above manner for the 12 stations revealed that the ob- 
served annual maximum events (regardless of type) at practically all stations 
are more frequent in the range of higher floods than is indicated by the com- 
ems combination frequency curves. This may be due to the fact that the _ : 
March 1936 and even the November 1927 non-hurricane floods are much more 
rare events than would be expected to occur in the 40-year period studied. _ 20 


Also, ‘the hurricane flood of August 1955 is considered to be much rarer than 
can be expected in the 40 years studied. On the other hand, it is possible that 
the log-normal technique employed, although successful in most regions of the 
United States, is not applicable in New England. Consideration was given to ae 
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e types, Dut Was concluded that other 1100d types were no 
— tinct or sufficiently independent. 
‘If the results of the study described herein were the only evidence available 
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the country, Py it is only reasonable to expect : a serious excess of large ¢ 
floods in some regions and a serious shortage of large floods in other re- aS : 
_ gions. _ As a matter of fact, there is a feeling among hydrologists familiar ' 
wit the New England region that the floods of 1936, 1955, and even 1927 are 
me very rare events, particularly in the light of historical data. It is be- 
lieved, therefore, that the application of the log-normal technique to the fre- — 
- quency analysis of hurricane floods and non-hurricane floods separately « con- 
Shoe, a reasonable method ad evaluating the future frequency of extreme 
floods in the New England area. The frequency curves derived by this 
method indicate that the 1955 floods were moderate at 9 of the 12 stations 7 
"studied, but that the flood of 18 August (Hurricane Diane) will be exceeded 4 
once in 200 to 300 years at the remaining three stations. 


pars 


— 
g the adequacy rong On 
Bee 
— 
— 


To). 


a 
= 
— 


LLOYD: 


